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I We present and develop a recursion scheme to construct joint eigenfunctions 

for the commuting analytic difference operators associated with the integrable N- 
particle systems of hyperbolic relativistic Calogero-Moser type. The scheme is based 
on kernel identities we obtained in previous work. In this first paper of a series we 
present the formal features of the scheme, show explicitly its arbitrary- N viability 
■ for the 'free' cases, and supply the analytic tools to prove the joint eigenfunction 

| properties in suitable holomorphy domains. 
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1 Introduction 

This paper is the first in a series of articles that are concerned with the explicit diag- 
onalization and Hilbert space transform theory for the relativistic generalization of the 
hyperbolic iV-particle Calogero-Moser system. As is well known, the nonrelativistic ver- 
sion is defined by the Hamiltonian 

fc2 N 

H 2 = --Y,dl ] +g{g-h) £ /i 2 /4sinh 2 (M^ - ^)/2), (1.1) 

j=l l<j<l<N 

where h is Planck's constant, g > a coupling constant with dimension [action], and fi > 
a parameter with dimension [position] _1 . There exist N — 1 additional independent 
PDOs Hk of order k, k = 1, 3, . . . , N, such that the PDOs form a commutative family. 
The simplest Hamiltonian is the momentum operator 

TV 

H x = -ih^2d x ., (1.2) 

i=i 

but the remaining Hamiltonians will not be specified here. (They can be found for example 
in the survey [QP83].) 

The arbitrary-g joint eigenfunctions of these PDOs were introduced and studied by 
Heckman and Opdam [HO 8 7] , and their Hilbert space transform properties were obtained 
by Opdam |Opd95| . (More precisely, these authors handle arbitrary root systems, whereas 
we restrict attention to Ajv-i-) For the case N = 2 the joint eigenfunction amounts to 
a specialization of the hypergeometric function 2 Fi, and the associated Hilbert space 
transform is a special Jacobi function transform, cf. Koornwinder's survey in |Koo84j . 

The relativistic generalization [RS86,Rui87] is given by the commuting analytic dif- 
ference operators (henceforth AAOs) 

Sk(x) = ^2 II f-( x m-x n ) JJexp(-z/i/3<9 ;C; ) Yl f+(x m -x n ), k = 1, ...,JV, (1.3) 



ic{i....,N} mei lei mei 

\i\=k n i l n £! 



where 



f ± (z) = (sinh(/i(z±z^)/2)/sinh(/iz/2)) 1/2 . (1.4) 

Here, > can be viewed as 1/mc, with m = 1 the particle mass and c the speed of 
light. In the nonrelativistic limit c — > oo these operators give rise to the above commuting 
PDOs. (See |Rui94j for a survey of the relativistic Calogero-Moser systems and their 
various limits.) 

Thus far, only for N = 2 the eigenfunctions and Hilbert space transform are well un- 
derstood. Indeed, they can be obtained by specializing results by the second-named au- 
thor pertaining to a 'relativistic' hypergeometric function generalizing 2 Fi [Rui9 9jlRui03II[ 
IRui03III] . This function is defined in terms of an integral whose integrand involves solely 
products of the hyperbolic gamma function from [Rui97j . (See also (vdB06, BRS07J for 
other perspectives on this function.) 
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In recent years, novel integral representations of the pertinent Ai-type (one-coupling) 
specializations of the latter BC\ (four-coupling) function have been obtained [Ruillj . 
These representations amount to Fourier transforms of products of hyperbolic gamma 
functions. For our purposes, the latter Fourier transform representations are of crucial 
importance. Indeed, as we shall show, they can be viewed as the result of the step 
from iV=ltoiV = 2ina recursive construction of the arbitrary- iV joint eigenfunctions 
of the A AOs S& ( 11.31) . The point is that the plane wave in the integrand can be viewed 
as the N = 1 eigenfunction, whereas the product of hyperbolic gamma functions serves 
as a kernel function, connecting the free N = 1 AAO exp(—i/3hd/dx) to the interacting 
N = 2 AAOs. 

In a recent joint paper [HRllj . we presented a comprehensive study of kernel functions 
for all of the Calogero-Moser and Toda systems of A N _i type. In particular, we obtained 
kernel functions connecting the hyperbolic AAOs for the iV-particle case to those for 
the (N — l)-particle case (see also [KNS09] ). For the case N = 2, the pertinent kernel 
functions amount to those occurring in |Ruill] , and this enables us to set up a recursion 
scheme for arbitrary N, as sketched in Section 2. 

The idea that such a recursive construction might be feasible is not new. It appears 
to date back to work by Gutzwiller [Gut81j . who used it to connect eigenfunctions for 
the periodic and nonperiodic (nonrelativistic) Toda systems. This formalism was then 
used for several other cases, in particular by Kharchev, Lebedev and Semenov-Tian- 
Shansky |KLS02] for the relativistic Toda systems, and by Gerasimov, Kharchev and 
Lebedev [GKL04] for the g = 1/2 specialization of the nonrelativistic hyperbolic Calogero- 
Moser system (cf. fll.ip ) and for the nonrelativistic Toda systems. We would like to express 
our indebtedness to this previous work, without which the scheme might seem to come 
out of the blue. 

In the later work following Gutzwiller's pioneering contribution, the representation 
theory of Whittaker modules and Yangians plays a pivotal role, and accordingly a consid- 
erable algebraic machinery occurs. However, as will become clear for the present case, the 
formal structure of the recursion scheme can be understood quite easily. Indeed, in our 
approach the main algebraic input consists solely of the kernel identities from our previous 
paper, and the same reasoning applies to the nonrelativistic hyperbolic Calogero-Moser 
Hamiltonians and to the nonperiodic Toda Hamiltonians, for which the pertinent kernel 
identities were also obtained in |HRllj . 

On the other hand, the simplicity of the construction in Section 2 hinges on glossing 
over a great many analytical difficulties. In fact, it is a major undertaking to show that the 
integrals yield meromorphic joint eigenfunctions that give rise to a unitary eigenfunction 
transform with the long list of expected symmetry properties. The snags at issue are 
already considerable for the first steps, and will become clear in due course. In this 
paper our focus is on a complete proof of the joint eigenfunction properties in suitable 
holomorphy domains, leaving various issues (including global meromorphy) open for now. 

We plan to address the analogous problems for the nonrelativistic hyperbolic Calogero- 
Moser and nonperiodic Toda systems in later papers. In particular, in the aforementioned 
work dealing with recursive eigenfunctions for these systems, the expected duality prop- 
erties are not shown and unitarity ('orthogonality and completeness') is left open. Also, 
the associated scattering theory needs to be studied, so as to confirm the long-standing 
conjecture that the particles exhibit soliton scattering (conservation of momenta and 
factorization). These features are quantum analogs of classical ones exhibited by the 
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action-angle transforms for these systems [Rui88j, and their relevance for the relation to 
the sine-Gordon quantum field theory has been discussed in |Rui01j . 

There are two length scales in the relativistic Hamiltonians (II. 3p . which we reparametrize 

as 

a + = 2tt/[jl, a_ = h/3. (1.5) 
Also, we trade the coupling g for a new parameter b with dimension [position], namely, 

b = (3g. (1.6) 

With these replacements in (11.31) in effect, the Hamiltonians with a + and a_ interchanged 
commute with the given ones, since the shifts change the arguments of the coefficients by 
a period. The resulting 2N commuting Hamiltonians can be rewritten as 

H k ,s{b\ x) = ] J fs,-{x m - x n )Y\_ew(-i a -sd Xl ) fs,+{x m - x n ), (1.7) 

IC{l,...,N}m£l lei mei 

\I\=k n n£I 



ss(z ± ib) \ 



where k — 1, . . . , N, S — +, — , and 

ss{z) 

Here and throughout the paper, we use the functions 

Ss(z) = smh.(nz/as), c$(z) = cosh(7r z/a$), eg(z) = exp(7iz/a$), 5 = +,—. (1.9) 
It is also convenient to use the parameters 

a = 27r/a + a_, a = (a+ + a_)/2, (1.10) 

a s = min(a + , a_), a; = max(a + , a_), a + ,a_>0. (1-H) 

To be sure, there are a great many different Hamiltonians commuting with the given 
Hamiltonians Hi + ...,H N+ . Indeed, one can replace the functions f-,±(z) in H k _(x) 
by arbitrary functions with period ia_. Of course, in that case the resulting Hamilto- 
nians H 1 , . . . ,H N will not pairwise commute any more. But the latter feature can be 
ensured by choosing any parameter b' that differs from b, and then all of the 2N Hamil- 
tonians do commute. Unless b' equals 2a — b, however, it it extremely unlikely that joint 
eigenfunctions of the Hamiltonians and H k exist. The choice b' = 2a — b is ex- 

ceptional, since it yields again the Hamiltonians H k _, as becomes clear by pushing the 
functions on the right of the shifts to the left, with arguments shifted accordingly. 

From the perspective of Hilbert space (which we do not address in this paper), it is 
crucial to restrict the parameters to the set 

n = {(a + , a_, b) E (0, oo) 3 | b < 2a}. (1.12) 

Clearly, whenever the coupling b is real, the 2N Hamiltonians ( 11. 71) are formally self- 
adjoint, viewed as operators on the Hilbert space L 2 (M JV , dx). To promote them to bona 
fide commuting self-adjoint operators, however, the restriction of b to the bounded in- 
terval (0, 2a) is already imperative for N = 2, since this key feature is violated for 
generic b > 2a, cf. [RuiOOj. 
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At this point we are in the position to add some further remarks about related litera- 
ture. First, there is Chalykh's paper |Cha02] . where Baker- Akhiezer type eigenfunctions 
of the above iV-particle Hamiltonians are introduced. With our conventions, these corre- 
spond to the special 6-choices b = ka + or b = fca_ with k integer. We intend to clarify the 
relation of the arbitrary-6 eigenfunctions furnished by the present method to Chalykh's 
eigenfunctions in later work. 

Secondly, there are several papers where so-called Harish-Chandra series solutions of 
the joint eigenvalue problem for the Hamiltonians are studied. A comprehensive study 
along these lines with extensive references is the recent paper by van Meer and Stok- 
man |VMS09j . Roughly speaking, in this setting one arrives at eigenfunctions that corre- 
spond to only one of the two modular parameters 

q + = exp(ma + /a_), g_ = exp(ina_/a + ), (1.13) 

on which our eigenfunctions depend in a symmetric way, in the sense that there is depen- 
dence on a single parameter q that must have modulus smaller than one. Accordingly, 
one only considers the above Hamiltonians H k j for one choice of 5. (For the BC\ case the 
relation between the latter type of eigenfunction and the modular-invariant relativistic 
hypergeo metric function has been clarified in [BRS07J.) 

We have occasion to use two further incarnations of the Hamiltonians Hk,s, viewed 
again as acting on analytic functions. These are obtained by similarity transformation 
with a weight function and a scattering function. The latter are defined in terms of 
the hyperbolic gamma function G(a + ,a_; z), whose salient features are summarized in 
Appendix A. 

First, we define the generalized Harish-Chandra function 

c(6; z) = G(z + ia- ib)/G(z + ia), (1.14) 
and its multivariate version 

C(b;x) = Yl c(b] Xj -x k ). (1.15) 
l<j<k<N 

(Here and below, we usually suppress the dependence on the parameters a + and a_; the 
dependence on b is often omitted as well.) Then the weight and scattering functions are 
defined by 

w{z) = l/c{z)c{-z), W{x) = l/C(x)C{-x), (1.16) 
u(z) = -c(z)/c(-z), U(x) = (-) N(N - 1 ^ 2 C(x)/C(-x). (1.17) 
Now we introduce 

A k , s (x) = W(x)- l / 2 H K5 (x)W(xf'\ (1.18) 

A k , s (x) = U{xy l l 2 H kt5 {x)U{xf/ 2 = C(x)- l A K5 (x)C(x), (1.19) 

where k = 1, . . . , N, and 5 = +, — . 

From the difference equations (1A.2D satisfied by the hyperbolic gamma function it 
follows that we have 

1 ileM-ia- S d Xm )W(x) = H - 2 f ( ^ X _ m ~li -y 



W(x) 1 A m A f fsJxm -x n - ia-s 

mdl mdl "it 
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Hence we obtain 



/c{i,..,iv}mej Wm XnJ lei 

\I\=k n 

Likewise, we deduce that the second similarity transformation yields 

a i \ TT s s (x m - x n - ib) s s (x m - x n + ib - ia- S ) T-f , ■ on 
A k ,s{x)= 2^ [I r — — l[exp(-ia_ s d Xl ). 

ic{i,...,N} meiMi Ss{Xm Xn) S6{Xm Xn ia ~ 5) lei 
\i\=k m > n 

(1.22) 

Thus the similarity-transformed AAOs act on the space of meromorphic functions. For 
parameters in II and x G K , the weight function W(x) is positive and the 'S'-matrix' 
U(x) has modulus one. Accordingly, the AAOs A ki $ and A k> s are then formally self- 
adjoint, viewed as operators on the Hilbert spaces L 2 (R N , W(x)dx) and L 2 {R N , dx), resp. 
(Once more, these features still hold for b real.) 
The scattering function satisfies 

U(2a-b;x) = U(b;x), (1.23) 

whereas the weight function is not invariant under the reflection b — > 2a — b. Therefore, 
the AAOs Ak t s are not invariant, whereas we have 

A kj5 (2a-b;x) = A k , s (b;x), k = l,...,N, 5 = +,-. (1.24) 

On the other hand, W(x) is symmetric (invariant under arbitrary permutations), whereas 
U(x) is not symmetric. (Indeed, w(z) is even, whereas u(—z) equals l/u(z).) As a 
consequence, the AAOs A k s are symmetric, whereas the A k ,s are not. (Note this can be 
read off from the restriction m > n in (11.221) .) More precisely, these different behaviors 
hold true for k < N — 1, since we have 

N 

H NtS = A NjS = A N ,s = J|exp(-za_ 5 a ci ). (1.25) 

i=i 

The 6-choices a + and a_ have a special status, inasmuch as they lead to constant 
coefficients in H k $ and A k> s (but not in A ki $). Indeed, we have 

H kj ±(a s ;x) = A k ,±{a s ;x) = 2J Il ex P(-* a =F^i)) k — l,...,N, 5 = +,-, (1.26) 

ic{i,...,N} lei 
\i\= k 

in accordance with no scattering taking place: 

U(a 5 ;x) = l, 6 = +,-. (1.27) 

(This follows from (11.141) and (11.171) by using (IA.2I) .) For these free cases, we shall show 
that the recursion scheme gives rise to the multivariate sine transform, and all of the 
expected properties can be readily checked. Even for these cases, however, some non- 
obvious identities emerge. This is because the kernel functions are already nontrivial for 
the free cases, and they are the building blocks of the eigenfunctions. 
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We proceed to sketch the results and the organization of the paper in more detail. 
In Section 2 we introduce the relevant kernel functions and their salient features, and 
present the recurrence scheme in a formal fashion (i.e., not worrying about convergence 
of integrals, etc.). As will be seen, the key algebraic ingredient for getting the eigenvalue 
structure expected from the explicit solution of the classical theory |Rui88] is given by the 
following recurrence for the elementary symmetric functions S% of M nonzero numbers 
ai, . . . , a,M '■ 

^fe M ^( a i) ■ ■ ■ i °m) = a M (*Sfc M l \ a i/ a Mi ■ ■ ■ i Qm-i/om) + S&^i ^(cli/o-m, • • • , a-M-i/ &m) J ■ 

(1.28) 

Here we have M > 1, k = 1, . . . , M, and 

fig f - 1) = 0, Sj> M - l) = l. (1.29) 

Section 3 is concerned with the free cases b = a±. For these cases the kernel functions 
reduce to hyperbolic functions, for which it is feasible to evaluate the relevant integrals 
explicitly. (A key auxiliary integral is relegated to Appendix C, cf. Lemma C.l.) Thus, 
the joint eigenfunctions can be obtained recursively, yielding the expected results. 

In Section 4 we focus on the analytic aspects of the first step of the scheme, allowing b in 
the strip Re b G (0, 2a). This step leads from the free one-particle plane wave eigenfunction 
to the interacting two-particle eigenfunction, and yields the relativistic conical function 
from |Ruill] (after removal of the center-of-mass factor). We reconsider some properties 
of this function, using arguments that do not involve the previous representations of the 
BC\ case (for which no multivariate recurrence is known). 

More specifically, we focus on holomorphy domains, the joint eigenvalue equations, 
and uniform decay bounds, with our reasoning (as laid down in Props. 4.1-4.4) serving as 
a template for the N > 2 case. In this special case, however, we can proceed much further 
than for N > 2. More precisely, we can easily obtain a larger holomorphy domain and 
corresponding bounds (cf. Props. 4.5 and 4.6), since contour deformations do not lead to 
significant complications. 

Section 5 is devoted to the step from N = 2 to iV = 3. This leads to novel difficulties, 
but the counterparts of Props. 4.1-4.4 can still be proved. To control contour deforma- 
tions, however, is already a quite arduous task for N = 3, and we cannot easily get the 
expected holomorphy features in this way. (In later papers we hope to clarify the global 
meromorphy character in both x and y for arbitrary N in other ways.) We do extend the 
holomorphy domain in the variable x (as detailed in Prop. 5.5), but the present method 
seems too hard to push through for arbitrary N. 

On the other hand, once our arguments yielding Props. 5.1-5.4 are well understood, 
the remaining difficulties for the inductive step treated in Section 6 are largely of a com- 
binatoric and algebraic nature. This relative simplicity hinges on the explicit evaluations 
of some key integrals, cf. Lemmas C.2 and C.3. It came as an unexpected bonus of the 
free case study in Section 3 that the integrals arising there (as encoded in Lemma C.l) 
suggested to aim for bounds involving the related integrals of Lemmas C.2 and C.3. In- 
deed, the latter furnish the tools to control the inductive step, which is encapsulated in 
Propositions 6.1-6.4. 

As we have mentioned already, this paper is the first in a series of articles. In Section 
7 we provide a brief outlook on future work. We discuss the main aspects of the joint 
eigenfunctions that we plan to investigate, and also mention some of the results we expect. 
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In Appendix A we review features of the hyperbolic gamma function we have occasion 
to use. In Appendix B we derive bounds on the G-ratio featuring in the kernel functions, 
and on the weight function building block w(z), cf. (11.14p — (11.16j) . 

As already mentioned, the explicit integrals needed to handle the free cases in Section 3 
can be exploited to explicitly evaluate two further integrals that are of crucial importance 
for the method we use to render the scheme rigorous. Indeed, the latter integrals enable 
us to derive in a recursive fashion uniform decay bounds on the joint eigenfunctions, which 
are needed to control the inductive step. Lemmas C.1-C.3 contain the statements and 
proofs of the pertinent integrals. 



2 Formal structure of the recursion scheme 

We begin this section by detailing the various kernel functions and identities. First, the 
special function 

^fj G ( x i ~Vk + ib/2) 

satisfies the kernel identities 

(A ktS (x)-A kJ {-y))S(x,y) = 0, k = l,...,N, 5 = +,-, (2.2) 



so that the functions 

and 

satisfy 



*{x,y) = [W{x)W{y)\ 1 ' 2 S{x,y), (2.3) 
K(x,y) = [C(x)C(-y)]- l S(x,y), (2.4) 



(H kiS (x)-H k , s (-y))*(x,y) = 0, k = l,...,N, 5 = +,-, (2.5) 

(A k , s (x) - A k>s (-y)))C(x, y) = 0, k = l,...,N, 5 = +,-. (2.6) 

(See |Rui0 6] for the proof of f!2 .2j) : the elliptic regime handled there is easily specialized 
to the hyperbolic one.) 

The kernel functions just defined connect the iV-particle AAOs to themselves. As we 
intend to show in a later paper, they have a rather surprising application to the study 
of the joint eigenfunctions produced by the scheme. The protagonists of the scheme, 
however, are kernel functions connecting the iV-particle AAOs to the (N — l)-particle 
AAOs, obtained in [HRllj . They arise from the previous ones by first multiplying by a 
suitable plane wave and then letting y n go to infinity. From now on, the dependence of 
the AAOs and kernel functions on N shall be made explicit wherever ambiguities might 
otherwise arise, and we also use a superscript fl to denote a kernel with first argument 
in C M and second one in C M ~ X . 

With these conventions in place, the kernel function 



s 



satisfies the key identities (cf. Corollary 2.3 in |HR11] ) 

A i k N s \x 1 ,...,x N )S t N {x,y) 

= {4 N s - 1 \-y 1 ,...,-y N „ 1 )+Al N sJ(-yi,...,-^^^ (2.8) 

where k = 1, . . . , N, 5 = +, — , and 

A^i-yu . . . , -y N ^) = 0, A^ (- Vl , . . . , -y N -x) = 1. (2.9) 

Using notation that will be clear from context, we now set 

¥ N (x,y) = [Wn^Wn-M^S^v), (2-10) 

4^,2/) = [C N (x)C N .. 1 (-y)]- 1 Si f (x,y). (2.11) 
The counterparts of are then (cf. also (1231)- (1231) ) 

fliy^^irCx.y) = (^ _1) (-y) + ^-M 1) (-y))^(^y), (2-12) 

4$W4r(*.v) = + 4 N rJ(-y))£ l N &yl (2.13) 

where fc = 1, . . . , N, S — +, — , and 

<-"(-!/) =<-"(-!/) -0, C(-,) = ^(-,) S L (2.14) 

We are now prepared to explain the 'calculational' crux of the recursion scheme. 
Assume we have a joint eigenfunction Jjv-i((^i, • • • , Xn-i), (yi, ■ ■ ■ , Vn-i)) of the AAOs 
A kS (x), with eigenvalues given by 

A{ k5^ l \ x ) J N~i{x,y) = S i k N ~ 1) (e s {2y l ), . . ., e s (2y N ^ 1 )) Jjv-i(x, y), k = 1, . . . , N - 1. 

(2-15) 

where S^\ai, . . . , cjm) denotes the elementary symmetric functions of the M numbers 
Oi, . . . , sm- Now consider the function Jn{x, y) with arguments x, y e C^, formally given 
by 

exp (iayN Ylj=i x j) f h 
Jn(x,v) = 7TT — TTi / d2W0v_i(2r)5^(x,2r)Jjv-i(2, (vi-Vn, ■ ■ ■ ,Vn-i-Vn))- 

(2.16) 

At this stage we do not address the convergence of the integral, and we also assume that 
we can take the x-dependent shifts in the AAOs under the integral sign. 

Acting with A k N s \x) on Jn, we pick up a factor esfikyx) upon shifting the AAO 
through the plane wave up front (recall (11.101) ). after which we act on the kernel function 
and use ( 12 .8p with y — )■ z. Using formal self-adjointness on L 2 (IR Ar_1 , W^-\{z)dz) of the 
two AAOs on the rhs, we now transfer their action to the factor Jjv-i, noting that the 
argument — z should then be replaced by z, since no complex conjugation occurs in ( I2.16p . 
As a consequence we can use our assumption (12.151) . the upshot being that we obtain 

A k N s( x ) J N(x,y) = e s (2ky N )[S i k N ~ 1) (e s (2(y 1 -y N ) ),..., e s (2 {y N _ x -y N ))) 

+ Sk-i 1] H2( Vl - y N )), . . . , e s (2(y N . 1 - y N )))} J N (x, y). (2.17) 



Invoking the symmetric function recurrence (ll.28p . the eigenvalue formula (I2.17P can now 
be rewritten as 

A${x)J N (x, y) = Sl N \e s (2 yi ), . . . , e 5 (2y N ))J N (x, y), k = 1, . . . , N. (2.18) 

Comparing this to our assumption (12.151) . we easily deduce that we have arrived at a 
recursive procedure to construct joint eigenf unctions. Indeed, we can start the recursion 
with the plane wave 

Ji(x,y) = exp(iaxy), (2-19) 

which obviously satisfies 

A$ k J 1 (x,y) = e±(2y)J 1 (x,y), (2.20) 
and then proceed inductively to obtain joint eigenfunctions for arbitrary N. 



3 The free cases b = a± 

The algebraic aspects of the procedure detailed in the previous section are easy to grasp 
and unassailable from a formal viewpoint. From the perspective of rigorous analysis, 
however, the scheme thus far has the advantage of theft over honest toil. Indeed, the first 
step in the recursion already leads to some delicate issues, as we shall see in the next 
section. 

On the other hand, the N — 1 starting point causes no difficulty Specifically, upon 
multiplication by (a + a_) -1 / 2 , the function (I2.19P yields the kernel of a unitary integral 
operator on L 2 (IR). Accordingly, the two analytic difference operators exp(— ia T d/dx) at 
issue can be promoted to commuting self-adjoint operators on L 2 {R,dx), defined as the 
pullbacks of the multiplication operators e±(2y) under Fourier transformation. 

Specializing to the case b = a+, we shall show in this section that the recurrence can 
be performed explicitly, yielding a multivariate version of the Fourier transform. (As will 
be seen shortly, the case b = a_ yields basically the same result.) 

We begin by writing down the specializations of the relevant functions. Setting b = a + 
in (12.71) . we deduce from the AAEs (1A.2|) satisfied by the hyperbolic gamma function that 
the kernel function S^ N reduces to 

N N-l 

Similarly, we find that the weight function Wjy reduces to 

W N (x)= Yl ^{Xj-Xkf, (3.2) 

l<j<k<N 

cf. (nn-dnsD. 

Turning to the explicit implementation of the recurrence, we start with the N = 2 
case, assuming x, y G M 2 until further notice. Substituting (I3.ip - (j3.2p into (I2.16P for 
iV = 2, we find 

r 2 i 

J 2 (x,y) = e ia ^ +X ^ / dzexp(iaz(y 1 -y 2 ))T[ 7 — -. (3.3) 
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Changing variable z — > a^z/rr and multiplying by 47r/a_, we see that for x\ 7^ x 2 and 
Hi 7^ Jj2 the resulting integral is equal to the right-hand side of ( 1C.1I) , with N = 1 and with 
t = 7rx/a_, pi = 2(yi — y2)/a + . Invoking Lemma C.l and making use of the identities 
cr(x) ■ y = x ■ o-~ l {y) and (— ) CT = (— ) fJ , we deduce that J 2 is given by 

J 2 (x,y) = - — r VV-^exp (iax ■ a(y)), (3.4) 

4s_(x 1 - x 2 )s + (y 1 - y 2 ) ^ 

where x,y el 2 and xi 7^ 0:2, ?/i 7^ ?/2- 

We claim that this structure persists in the general case, in the sense that 

M*,v) = II 4s TZ _ r ,T" 77. _ m) J2 (")' Jex P i iax ■ °(v))> ( 3 - 5 ) 

where x,y E M. N , and Xj 7^ ?/-,• 7^ for j 7^ A;. To verify this claim, we proceed by 
induction on N. Since the case N = 2 has just been shown, we assume (I3.5P for N > 2 
and prove its validity for N — > N + 1. 
From fl2TT5|) and (EH]) (EI2D , we infer 

exp(ia?/Ar + i(xi H Kxat+i)) t-i- -ia_ 



■W*>v) = oN( N+ i) A n 11 



2 N(N+1) N \ 11 s (y y) 

l<j<k<N +Kyj y/c7 



x Yl / rfzexp ( iQ! Z^^(^(i) -Vn+i) ) V — v ( 3 - 6 ) 

Now, we fix cr G S*^ and consider the corresponding integral in (13.61) . Changing variables 
zu — > a_Zk/iT, k = 1, . . . , N, and multiplying by (7i/a-) N , we arrive at the right-hand side 
of ( IC.lj) for t = 7rx/a_ and = 2{y a ^ — y^ + i)/a + , j = l,...,N. Hence, Lemma C.l 
yields 



J N +i(x,y) = — TT - — % -%— f r 

x EH ff E (-rexp^r^)-^)). (3.7) 

Here, we have identified a G SV with the element in Sjv+i that acts by a on the first N 
coordinates and leaves the last one invariant. Upon replacing r(s) ■ a(y) by x ■ (r" 1 ^)^/), 
substituting r — > ar^ 1 , and then using (— ) fTT = (— )°"(— ) T , we reduce the double sum- 
mation to 

N\ J2 (-) T exp(iax-T(y)). (3.8) 

If we now take r — > a, then we arrive at (13. 5p for iV — > N + 1. 
We now summarize and slightly extend our finding. 

Theorem 3.1. Assume b is equal to a + or to a_. Then we have 

J N {a s ;x,y)= TT V (-) CT exp (iax ■ a{y)), 5 = +,-, 

(3.9) 

where x, y G WL N , and Xj 7^ Xk, yj 7^ yu for j 7^ k. 



11 



Proof. We have just shown this equality for 5 = +. We need only interchange a+ and a_ 
throughout to obtain (13.9p for 5 = — . □ 

Using the G-AAEs flA.2j) . one easily checks that the definitions fll.14p - fll.15p entail 

C(a s ;x) = 1/ Y[ Hs-s&j - x k ). (3.10) 

l<j<k<N 

Thus it follows that the similarity transformations fll.l8l) - fll.l9p give rise to joint eigen- 
functions of the Hamiltonians H^ g^ag] x) and Ak,s{o>8', x) that are proportional to the 
kernel of the multivariate sine transform, 

^(x,y)= exp (iax ■ a(y)), (3.11) 

in accord with their free character for b = a±, cf. fll.26p . 



4 The step from N = 1 to N = 2 

In this section we consider the N = 2 case for 6-values in the strip 

S a = {b G C | Re 6 G (0,2a)}. (4.1) 
Fixing i/6l 2 until further notice, we begin by studying the integrand 

2 G(x — z — ib/2) 

h{b;x,y,z) = S s 2 (b; x, z)J 1 (z,y 1 - y 2 ) = exp(iaz(y! - y 2 )) \\ Q ^ _ - - , (4.2) 

arising in the first step J\ — > J 2 , cf . (I2.16P . In the z-plane it has upward / downward double 
sequences of G-poles at (cf. fIXTTTTt — <IX7T2j> ) 

z — Xj — ib/2 + ia + zm, z = Xj + ib/2 — ia — zm, j — 1, 2, k,l e N. (4.3) 

Letting first x G 1R 2 , it is clear that the integration contour M. stays below/above these 
upward/downward sequences. From Lemma B.l we see that I 2 decays exponentially for 
z —> ±00, so the function 

J 2 (6; x, y) = exp(iay 2 (xi + x 2 )) \ dzl 2 (b; x,y,z), be S a , x, y G M 2 , (4.4) 

is well defined. 

We proceed to obtain holomorphy properties of J 2 in the variables b and x, using 
arguments that allow a recursive generalization to arbitrary N. (For the N = 2 case, 
it is easy to obtain more information, as will transpire after Prop. 4.4.) First, from 
the fixed contour representation f)4.4p and pole locations (14.31) . it is immediate that J 2 
is holomorphic in b and Xj for b G S a and |Ima;j| < a — Re 6/2, j = 1,2. Secondly, 
Lemma B.l entails 

^2(6; x, z) = 0(exp(— aReb |Rez|)), |Rez| — > 00, (4.5) 
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where the implied constant is uniform for (6, x, Im z) vaying over compact subsets of 
S a x C 2 x R. Thus we can freely shift the contour R up and down as long as we do not 
meet any poles. 

In particular, for x G R 2 we may replace the contour R in (I4.4p by a contour R — za + 
z'6/2 + ze, e G (0,a — Re 6/2). Then we can read off that J 2 extends holomorphically to 
Imxj G (e, —2a + Re 6 + e), j = 1,2. Clearly, this argument can be iterated. Likewise, 
we can move the contour up step by step to allow arbitrary positive Imxj with similar 
restrictions. More precisely, introducing the domain 

D 2 = {(b,x) ES a xC 2 \ |Im(xi - x 2 )\ < 2a -Re b}, (4.6) 

we have the following result. 

Proposition 4.1. Let y G R 2 . Then the function J 2 (b;x,y) is holomorphic in D 2 . 
Proof. We fix (6, x) G D 2 and set 

r] = Im(xi +x 2 )/2. (4.7) 

Then we have 

Imxj + a - Re 6/2 > 77, Imxj — a + Re 6/2 < 77, j = 1,2, (4.8) 

so the contour R + ir\ stays below/above the upward/downward G-pole sequences (14. 3p . 
Therefore, we can arrive at these x-values by the above iterative contour shift procedure 
without passing any poles. Hence we deduce holomorphy near (b,x), and thus in D 2 . □ 

The analyticity properties obtained thus far already suffice to show that J 2 is a joint 
eigenfunction of the N = 2 AAOs with the expected eigenvalues. More precisely, we fix 
attention on the subset 

D 2 . s = {(b,x) G D 2 I (6, (xi - ir)i,x 2 - irj 2 )) G D 2 ,\/(r)i,r] 2 ) G [0,a ; ] 2 }, Re 6 G (0,a s ), 

(4.9) 

for which the argument shifts stay in the holomorphy domain D 2 . Here we have used 
the notation f ll.lip . noting that the restriction of Re 6 to (0, a s ) is necessary for D 2 s to 
be non-empty. (Once we obtain further analyticity properties in 6 and x, the eigenvalue 
equations can be analytically continued to 6 G S a , cf. below.) 

Proposition 4.2. Let y G R 2 . For all (b,x) G D 2)S , we have the joint eigenfunction 
property 

A^(x)J 2 (x,y) = S^(e s (2 yi ),es(2y 2 ))J 2 (x,y), k = 1,2, 5 = +,-. (4.10) 

Proof. Thanks to analyticity of J 2 in D 2 and the restriction to D 2jS , the shifts of Xj 
by —ia_s are well defined. Also, by analyticity in the strip Re 6 G (0, a s ), it suffices 
to prove (14.101) for 6 G (0, a s ). Then the distance a s + a\ — 6 between the upward and 
downward pole sequences in (14.31) is larger than a\ by e = a s — 6 > 0. Moreover, by 
analyticity in x, we need only show (14.101) for |ImXj| < e/2, j = 1,2 (say). 

Letting x vary over the region just detailed, we now proceed in a number of steps to 
rewrite the lhs of (14.101) . First, we note that for this x- region the ^-contour R in (14.41) is 
not only permissible, but can be moved down to 

C^=R-iai/2, (4.11) 
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without meeting any poles. Second, after shifting the AAO through the plane wave up 
front (picking up a factor es(2ky 2 )), we are allowed to act with it under the integral sign, 
since the poles stay away from C_ while they are shifted down by a distance as- Thus 
the lhs equals 

e s (2ky 2 )e iay ^ + ^ [ dzA®(x)S i 2 (x 1 z)J 1 {z,y 1 -y 2 ). (4.12) 

JC- 

Third, we use the kernel identity (12. 8p to trade the action of Af\(x) on S 2 (x, z) for that 
of exp(iasd/dz) for k = 2 and of exp(ia_sd/dz) + 1 for k — 1. Then the integral in (I4.12p 
can be written as 

dzSl(x,z + ia- 5 )J 1 (z,y 1 - y 2 ) + (2 - k) / dzl 2 (x,y,z). (4-13) 
c_ JC- 

Fourth, we change variables z — > z — ia_s in the first integral of (14.131) . which yields 



dzSl(x,z)Ji(z — ia-s,yi — y2) + (2 — k) / dzl 2 (x,y } z). (4.14) 

C-+ia_ s JC- 

Fifth, we use the eigenvalue equation (recall ( I2.19p ) 

J x {z - ia_j, Vi - Jfe) = e 5 (2(?/i - y 2 ))J\{z, y x - y 2 ), (4.15) 

and pull out the eigenvalue from the integration. Sixth and last, we shift the contours 
C_ + ia_s and C_ (for k — 1) back to R without meeting poles, obtaining ( 14.10]) . Thus 
the proof is complete. □ 

Next, let us take z — > z + (xi + x 2 )/2 in the integral ( 14. 4ft . This yields 
J 2 (&; (xi,x 2 ), (yi,V2)) = exp(za(xi + ar 2 ) (2/1 + y 2 )/2)J(b; x 1 - x 2 ,y 1 - y 2 ), (4.16) 
where (recall the reflection equation (IA.6|) ) 



J(b; u, v) 



( dze iazv I [ G(5 lZ + 5 2 u/2-ib/2), u,v e R, 6 G S* a . (4.17) 



From this it is clear that J is even in u and v. Furthermore, it follows that J 2 is invariant 
under the interchanges X\ <H- x 2 and y\ ^ y 2 . Note that the first symmetry property is 
already plain from the defining formula (I4.4p . whereas the second one is not immediate 
from flOp . 

Comparing (I4.17P to Eq. (3.51) in [Ruillj . we see that J is related to the relativistic 
conical function studied in (Ruillj via 

J(b; x, y) = y/a + a_ 'G{ia-2ib)n{b;x,y) JJ G(5y-ia + ib). (4.18) 

s=+ - 

From this a rather complete picture of J and J 2 can be deduced. In order to achieve 
analytic control on the second step in the scheme, however, we need a uniform bound on 
the 2-particle eigenfunction J 2 , which exhibits its exponential decay for |Re (x\ — x 2 )\ — > 
00. This bound is not immediate from the results obtained in [Ruillj . 
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Proposition 4.3. Let y G M 2 . For any (b,x) G D 2 we have 

\J 2 (b;x,y)\ < C(b,\lm(x 1 -x 2 )\)exp(-alm(x 1 +x 2 )(y 1 + y 2 )/2) 

x . ^\~ X2) V V 7 = «Refe/2, (4.19) 
smn(7 Ke [Xi — x 2 )) 

where C is continuous on S a x [0, 2a — Re b). 

Proof. In view of ( 14. 16ft . this amounts to a bound on J(b;u,v). Clearly, the representa- 
tion (I4.17P can be analytically continued to the strip |Imw| < 2a — Re b, and for u in this 
strip we have 

\J(b;u,v)\< dz\r{b;z + u/2)r{b;z-u/2)\. (4.20) 
Jr 

From Lemma B.l we now deduce 

is cosh(7(2 - <5Re u/2)) 

where the C r -product is bounded above by a function C(b, \ lmu\) that is continuous on 
S a x [0, 2a — Re b). The remaining integral 

dz (4.22) 



/ 

JR 



cosh(7 (z — Re u/2)) cosh(7 (z + Rett/2)) ' 

can be easily evaluated by a residue calculation. (Alternatively, it follows from the N = 1 
case of Lemma C.I by letting p converge to zero.) Hence we arrive at ( I4.19p . □ 

Thus far we have assumed that y is real. As a consequence, the bound (14.201) is 
independent of v — y\ — y 2 - But when we combine (I4.16p - (l4.17p with (14.2 lj) . it becomes 
obvious that J 2 has an analytic continuation to all y G C 2 satisfying |Im (yi — y 2 )\ < Re b. 
Moreover, we can readily estimate J 2 in the larger holomorphy domain 

V 2 = {(b,x,y) G S a x C 2 x C 2 | (b,x) G D 2 , |Im (y\ — y 2 )\ <Re6}. (4.23) 

Proposition 4.4. For any (b,x,y) G T> 2 with Im (yi — y 2 ) ^ we have 

\J 2 (b;x,y)\ < C(b,\lm(x 1 -x 2 )\)exp(-alm[(x 1 +x 2 )(y 1 +y 2 )]/2) 

x sinh(alm ( Vl - y 2 )Re (xi - x 2 )/2) ^ 
sin(7rlm (y 1 — y 2 )/Re b) sinh(7 Re (xi — x 2 )) ' 

where 7 = aReb /2 and C is continuous on S a x [0, 2a — Re b). 

Proof. Once more, this is really a bound on J(b;u,v), cf. (14. 16j) (14. 17p . In this case we 
deduce from Lemma B.l that we have 

\J(b;u,v)\< [ dzeM-azlmv) J] ^^"^L . (^) 
Jr s ±±_cosh{7(z-5Reu/2)) 

so now we obtain 

f exp( — c^zlm v 1 

\J(b;u,v)\<C(b,\lmu\) / dz — — ; |Im^| < Reft, 

J R cosh(7 [z — Re u/2)) cosh(7 [z + Re u/2)) 

(4.26) 

where C is continuous on S a x [0, 2a — Re b). The integral can be evaluated by a residue 
calculation, or follows directly from Lemma C.I with N — 1. From this result we obtain 
the bound <^M). □ 
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Clearly, we can view Prop. 4.3 as a corollary of Prop. 4.4, obtained by letting Im (y± — 
y 2 ) converge to 0. We have distinguished the two settings for expository reasons. Indeed, 
we need only invoke Prop. 4.3 in the next section to obtain the N = 3 counterparts of 
Propositions 4.1 and 4.2, dealing with the joint eigenfunction properties for real y G M 3 . 
Then we obtain the N = 3 analogs of Propositions 4.3 and 4.4, which suffice for the 
inductive step taken in Section 6. 

For the cases N = 2 and iV = 3, we can readily derive further significant information, 
which we were unable to generalize inductively thus far. To start with, in the present 
N = 2 case we can readily obtain further analyticity properties of J 2 , and we continue by 
doing so. 

To this end, we fix y G C 2 with |Im (yi —7/2) | < Re 6. A moment's thought reveals that 
suitable contour deformations imply that J 2 has a holomorphic continuation to all (b, x) G 
S a x C 2 such that x\ — x 2 + ib and x 2 — X\ + ib stay away from the 'cut' i[2a, 00). Indeed, 
we need only indent the contour so as to ensure that the upward and downward pole 
sequences (14. 3 p remain above and below the contour, resp. Furthermore, contour pinching 
for ±Re (xi — x 2 ) = Imb can be avoided by requiring |Im (x± — x 2 ) \ < 2a — Re b. We now 
summarize this analysis. 

Proposition 4.5. The function J 2 (b;x,y) is holomorphic in (b,x,y) on the extended 
domain 

T>2 = {(b,x,y) G S a xC 2 x C 2 | ±(x x -x 2 )+ib <£ i[2a,oo), \1m(y 1 -y 2 )\ < Re b}. (4.27) 

As announced, the AAEs (I4.10p can now be analytically continued to the subset 

V\ s = {(b,x,y) G VI I {b,{x 1 -ir ]1 ,x 2 -ir] 2 ),y) G V e 2 ,\f(VuV2) e [0,a,] 2 }, (4.28) 

for which the argument shifts stay in the holomorphy domain T> 2 . In particular, we have 

{(b,x,y) G (0,2a) x C 2 x R 2 | Rexi ^ Rex 2 } C V e 2 s . (4.29) 

Shifting the contour beyond the poles of the integrand, it is possible to show global 
meromorphy in x, but to detail this is beyond the scope of the present paper. We do 
want to round off this section, however, by deriving a generalization of Prop. 4.4 to the 
extended holomorphy domain Pf. By now it will be clear that this amounts to a further 
bound on J(6; u, v) on its holomorphy domain 

Vj ee {(b,u,v) G S a x C x C I ±u + ib <£ i[2a,oo), |Imv| < Re 6}, (4.30) 

cf. (insHnzi). 

Proposition 4.6. For all (b,u,v) G T>j we have 

1 TV 7 \ 1 \ sinh(aImwReM/2) ... 

J{b;u,v) <C{b,u,v) , ' 1 (4.31) 

sm(7rlmf/Keo) smh(7HewJ 

Here, C(b,u,v) is a continuous function on T>j, which satisfies 

C(b,u,v) = 0(1), |Rew|^oo, (4.32) 

where the implied constant can be chosen uniformly for (6, Imu, Rev) varying over compact 
subsets o/S„xlxI 
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Proof. From holomorphy of J in T>j it is plain that a bound of the form ( I4.3ip is valid 
for |Rett| < R, with R > arbitrary and C(b,u,v) continuous on T>j. Fixing K > and 
choosing 

R=\lmb\+2a + K, (4.33) 

(say), it remains to obtain an estimate (14.311) for b G S a and u£C satisfying |Rett| > R 
and |Imtt| < K, with C(b,u,v) replaced by a function that only depends continuously on 
b, K and Rev. 

By evenness, we need only consider the case Rett > R. Then we can represent J as 



J(b;u,v) = J dz exp(iazv)r(b; z)r(b; z — u), (4.34) 

where C is a piecewise linear contour {s + ifxi 3 ) I s G E}, defined by setting 

/ K (s) = 0, s G (-oo,0] U [Rew + i^oo), (4.35) 

Ik^s) = Imw, s e [Rett — |Im6|, Rett + |Im6|], (4.36) 

and letting Jk{s) be a linear function so that C connects the points 0, u— \lmb\ and u + 
|Im6|, Rett + R in the z-plane. Since we require |Imtt| < K, we have |/^(s)| < 1 on C. 
On the slanting parts of the contour, the increase in length for a given increase in the 
curve parameter s is therefore bounded above by \/2 (compared to the flat pieces). As a 
consequence, we deduce a majorization 

\J(b;u,v)\ < V2exp(aK\Rev |) / ds exp(— aslmt>)|r(6; s + ifK(s))r(b; s + — u)\. 

Jr 

(4.37) 

Next, invoking the r-bound ( lB.6l) -( !B~7l) . we readily infer that the product 



C r (b,z)C r (b,z-u), zeC, (4.38) 

can be bounded above by a constant depending only on b and K. Thus the integral 
reduces to the one in (I4.26p . and so the assertions follow. □ 



5 The step from N = 2 to N = 3 

In this section we focus on the analytic aspects arising in the step ^ = 2—7-^ = 3, 
following the flow chart of Section 4 as far as possible. We recall the pertinent integrand 
is given by (cf. (12 . 161) ) 

hip; x, y, z) = W 2 (b; z)S\{b] x, z)J 2 (b; z, (y 1 - y 3 , y 2 - t/ 3 ))- (5.1) 

From Lemma B.2 we can now obtain bounds on the 2-particle weight function defined 
by (j!.14p — (11.16p . Likewise, Lemma B.l supplies bounds on the G-ratios in <S|, and 
Prop. 4.3 on J 2 . 

More specifically, letting y G M. 3 and Im {z\ — z 2 ) = until further notice, we obtain 
from (1B.6P 

l#i^)|<nf[ Crib^Zk-Xj) Reb/2 (5 . 2) 

fi fe = i cosh( 7 Re (z k - Xj)) 
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provided (6, Zk — Xj) G V r for j = 1, 2, 3, and k — 1, 2; from (IB . 13|) we have 

|W 2 (6;z)| < sinh 2 ( 7 (zx - z 2 )), Im^ = Imz 2 , (5.3) 

whereas (I4.19P yields 

1^2(6; z, (3/1 — j/3, 3/2 -2/3)) I < C(6)exp(-aIm(zi + z 2 )(yi + y2-2y 3 )/2) 

x . f 1 , — — — rr, Im^i = Im^. (5.4) 
smh(7(zi -2; 2 )J 

It is clear from these bounds that I 3 has exponential decay for |Re^i| — > 00 and for 
|Rez 2 | — > 00, uniformly for Im Zj varying over bounded intervals. 

Requiring at first x G R 3 , the integration contour R in the z^-plane stays away from 
the G-poles at (cf. flATT]) - fr02|) ) 

Zk = Xj — ib/2 + ia + zi m , z^ = Xj+ib/2 — ia — zi m , k = 1,2, j— 1,2,3, !,m£N, (5.5) 

so that 

e iay s (x 1 +x 2 +x 3 ) r 

J 3 {b; x, y) = / dzl 3 {b; x, y, z), (6, x, y)eS a xR 3 x R 3 , (5.6) 

2 JR2 

is well defined. 

It follows from the above that J 3 (as initially defined by f)5.6p ) extends to a function 
that is holomorphic in b and Xj for b G S a and |ImXj| < a — Re 6/2, j = 1, 2, 3. Further- 
more, thanks to the above bounds, we can shift the two contours R simultaneously up 
and down as long as we do not meet the poles (I5.5p . Using the same arguments as for the 
case N = 2, we can now extend the holomorphy domain step by step. In order to detail 
the result of this iterative procedure, we introduce 

D 3 = {(b.x) G S a x C 3 I max llm (x,- - x k )\ < 2a - Reb}. (5.7) 

l<j<fc<3 

Also, given x G C 3 , we set 

(f)(x) = lm(xj 1 +Xj 3 )/2, (5.8) 
where the indices ji,j 3 are such that 

Imx i3 < lmx h < lmx h , {ji,j 2 ,J 3 } = {1,2,3}. (5.9) 

Then we have the following counterpart of Prop. 4.1. 

Proposition 5.1. Let y G R 3 . Then the function J 3 (b;x,y) is holomorphic in D 3 . 

Proof. Fixing (b,x) G D 3 , we clearly have 

Imxj + a - Reb/2 > </>(x), Imxj - a + Reb/2 < (p(x), j = 1,2,3. (5.10) 

Therefore, the contour M + i(p(x) stays below/above the upward/downward G-pole se- 
quences. Hence we can continue J 3 to any such x-value by simultaneous contour shifts, 
without passing any of the poles (I5.5P in the iteration. More specifically, we arrive at the 
associated representation 

e iay3(xi+X2+xs) r 

J 3 (b;x,y) = / dzl 3 (b;x,y,z), (b,x) e D 3 , y G R 3 , (5.11) 

Z J (R+iij)(x)) 2 

from which holomorphy in D 3 can be read off. □ 
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Just as for the first step, these analyticity features are sufficient to show that J 3 is a 
joint eigenfunction of the six AAOs at hand. Accordingly, we introduce the subdomain 
of D 3 given by 

D 3)S = {(b, x) E D 3 I (b,x-irj) E D 3 ,V V e [0,a,] 3 , RefeG (0,a.)}. (5.12) 

On the other hand, we now have more than one contour, and in the proof of the follow- 
ing counterpart of Prop. 4.2 we have occasion to shift one of them. Thus, to show that this 
causes no problem with the contour tails, we can no longer appeal to the bounds (15.31) - 
(15. 4p . which require Im^ = Im2 2 . 

However, this snag is readily obviated, as follows. First, we note that the bounds ( IB. 141) - 
(1B. 1 5[) imply that the rate of divergence of W 2 (b; z) as |Re {z\ — z 2 )\ — » 00 does not depend 
on Im [z\ — z 2 ). Second, we may invoke Prop. 4.3 as long as |Im (z\ — z 2 )\ < 2a — Re b, and 
on closed subintervals the decay rate in the bound on the ^-factor in the integrand ^3 (15. ip 
as I Re (zi — z 2 )\ — > 00 again does not depend on Im (zi — z 2 ). In view of the bound (15.21) 
on the kernel function, this suffices for the contour shifts to be legitimate. 

Proposition 5.2. Fixing y G M 3 and letting (b,x) e D 3s , we have the eigenvalue equa- 
tions 

Af}(x)J 3 (x,y) = S^(es(2y 1 ) J es(2y 2 ),e s (2y 3 ))J 3 (x,y) } k = 1,2,3, 5 = +,-. (5.13) 

Proof. We follow the same line of reasoning as in the proof of Prop. 4.2. Thus we infer 
again that we need only prove (I5.13P for b G (0,a s ) and |ImXj| < e/2, j = 1,2,3, with 
e = a s — b. 

Letting x vary over this region, we now adapt the six steps in the proof of Prop. 4.2. 
First, the z^-contours ]R in (15.61) are permissible and can be moved down to C_ (14. lip 
without meeting poles. Second, shifting the AAO through the plane wave up front (yield- 
ing a factor es(2ky 3 )), we are allowed to act with it under the integral sign, so that the 
lhs of (I5.13P becomes 

e iay 3 (x 1 +x 2 +x 3 ) f f „ 

e 5 (2ky 3 ) d Zl dz 2 W 2 (z)Af}(x)Si(x, z)J 2 (z, ( yi -y 3 ,y 2 -y 3 )). (5.14) 

Third, we use (12. 8p to replace A ( k ' s (x) by the sum of A K k ' 5 (—z) and A k _ 15 (—z). Then it 
remains to show that the integrals 

/ dz x [ dz 2 W 2 (z)J 2 (z,( yi -y 3 ,y 2 -y 3 ))A^(-z)Sl(x,z), (5.15) 

JC- JC- 

are equal to 

Sf ] {e s {2(yi - y 3 )),e 5 (2(y 2 - y 3 ))) 2exp(-iay 3 (x 1 + x 2 + x 3 ))J 3 (x,y), (5.16) 
for / = 1 and I = 2. 

Fourth, we change variables depending on the two cases at hand. For / = 2 the AAO- 
action yields S^(x, (z\ + ia_$, z 2 + ia_s)) an d we change variables z^ — > 2% — ia_$, k = 1,2, 
yielding the integral 

/ dz\ I dz 2 W 2 (z)J 2 ((z 1 -ia_s,z 2 -ia- 5 ),(y 1 -y 3 ,y 2 -y 3 ))S 3 i (x,z). (5.17) 

J C-+ia_s JC-+ia_s 
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By virtue of Prop. 4.2 (with k = 2), this is equal to 



eg (2 {y x -y 3 ) +2 (y 2 -y 3 )) / dz 1 / dz 2 W 2 (z)J 2 ((z h z 2 ),(yi-y 3 ,y 2 -y 3 ))Sl(x,z), 

JC-+ia_s JC-+ia_g 

(5-18) 

which amounts to the fifth step for this case. Shifting the contours back to R (the sixth 
step), we see this equals (15 . 1 6[) for the case I = 2. 

It remains to show equality of (I5.15P and (I5.16P for 1 = 1. In this case the integral (I5.15P 
is given by 

dz l j dz 2 W 2 (z)J 2 (z,(y 1 -y 3 ,y 2 -y 3 )) 

C- JC- 

x (^—7 — Z2 exp(ia- S d zl ) + ^— ^ — Zl + ^ exp(ia_ 5 d Z2 )]sl(x, z), (5.19) 
V s s (zi - z 2 ) s s (z 2 - Zi) ) 



cf. f ll.2ip . Note that the pole at Z\ = z 2 due to the denominator s$ is matched by a double 
zero due to W 2 (z) = w(z\ — z 2 ), cf. ( Il.isp . 

Changing variables, we now rewrite (15. 19[) as 

, f , s s {zt - ia^s - z 2 + ib) . 
dzi / dz 2 -. ; r — W 2 (zi - ids, z 2 ) 

C-+ia_ s JC- Wl - * a -S ~ Z 2 ) 

x Sl(x, z) J 2 {{z 1 - ia^s, z 2 ), (y x - y 3 , y 2 - y 3 )) 

f , f , s 5 (z 2 - ia^s ~ zt + ib) . 

+ / dzt dz 2 ; r — W 2 (z 1 ,z 2 - ia- S ) 

Jc^ JC-+ia^ s s s (z 2 - ias - Zl) 

x Sl(x, z)J 2 ((zx, z 2 - ia_ 5 ), (y x - y 3 , y 2 - y 3 )). (5.20) 
From the AAEs (ll.20p satisfied by the IV-function we deduce that this equals 

dZl f dZ2 s ^-^-^ w 2 (z) 

C-+ia_ s JC- s S[Zl — Z 2 ) 

x S\{x, z)J 2 ({z 1 - ia_ s , z 2 ), (y x - y 3 , y 2 - y 3 )) 
+ / dzi dz 2 r — W 2 (z) 

JC- JC-+ia_s S S (Z 2 - ZI ) 

x Sl(x, z)J 2 ((z u z 2 - ia_ s ), (yi - y 3 , 2/2 - y 3 ))- (5.21) 

We claim that when we shift the contour C_ in these two integrals up by as, then no 
poles are met. (In view of the paragraph preceding (15.131) . we retain exponential decay on 
the contour tails, so the shift causes no problems at the tail ends.) Taking this claim for 
granted, the integrands are regular for zi, z 2 G C_ + ia_s, and since the double contour 
in the first integral is now equal to the one in the second integral, we are entitled to use 
the eigenvalue equation (I4.10p with k = 1 (the fifth step). The sixth step is the same as 
for the case / = 2 and now yields ( 15.161) for / = 1. 

To complete the proof, it remains to verify the claim. Clearly, we stay in the holo- 
morphy domains of 1S3 and J 2 while shifting the relevant contour. Setting z = z\ — z 2 
and letting first a^s = d.s, we meet the simple poles of the factors 1/G(±z + ia — ib) at 
±z = ib, but they are matched by zeros of s$(±z — ib). Likewise, the simple poles of ss{z) 
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at z = and at ±z = ia s (the latter being met when a s = ai) are matched by zeros of 
G(±z + ia). Thus the first integrand is regular for z\ G C_ + ia s and Im [z\ — z 2 ) G [0, a s ] 
and the second one for z 2 G C_ + ia s and Im (z 2 — Z\) G [0, a s \. 

Finally, consider the case a_$ = ai. Then we meet simple poles of 1/G(±z + ia — ib) 
at ±z = ib + ika s for fca s < a\ — 6, but they are again matched by zeros of ss(±z — ib). 
The poles of sg(z) at ±z = ika s for /ca s < ai are matched by zeros of G(±z + ia). Hence 
the first integrand is regular for z\ G C_ + z'a; and Im (zi — z 2 ) G [0, a/] and the second 
one for z 2 G C_ + z'a; and Im (z 2 — Zi) G [0, ai\. Thus our claim is proved. □ 

We proceed to obtain the analog of Prop. 4.3. Fixing (b, x) G D%, we recall the 
definition ( 15. 8 j) of 0(x), and define in addition (cf. (15. 9p ) 



d(x) = Im (xjj — Xj 3 ). (5.22) 
Proposition 5.3. Let?/ G M 3 . For any (6, x) G -D3 we have 



\J 3 (b;x,y)\ < C(b,d(x))exp(-a[(j){x){y 1 + y 2 ) +y 3 Imx j2 ]) 

x TT Re (x m - x n ) (5.23) 

. . , sinh(7 Re (x m - x n )) ' 

l<m<n<3 

where C is continuous on S a x [0, 2a — Re b). 
Proof. From ( 15. lift we obtain a majorization 

ii m ^ exp(-ay 3 Im (xi + x 2 + x 3 )) /" jiw^cS^ u/ / vu 

h\?,V)\< o / ^|W2(z)<S|(x,z)J 2 (^, (2/1—2/3, 2/2—2/3))!. 

«/(IR+M/>(:r)) 2 



(5.24) 



2 

Combining Prop. 4.3 with the bounds ( 15. 2 p and (15. 3p . we infer 

\J 3 (x,y)\ < C(b) exp(-a[y 3 Imx i2 + (y x + y 2 )<p(x)}) 

x / ds( Sl - S2 ) sinh (7(gl - 82)) n n Cr{ \?t " tt ' ( 5 - 25 ) 

J M2 J= ifc=i cosh( 7 (s fc - Rexj)) 

with C(6) continuous on S a . Now we have (cf. (15. 8p . (I5.22j) and (15. 7p ) 

|Im(i0(x) -Xj)\ < d{x)/2 < a -Re 6/2, j = 1,2,3, x G D3. (5.26) 

Recalling Lemma B.l, it follows that the C r -product is bounded above by a function 
C(b, d(x)) that is continuous on S a x [0, 2a — Re 6). The remaining integral 

ds (^72)smh( 7 ( 3l - 32 )) (52?) 
n i= i Ilfc=i cosh( 7 (s fc - Rex^) 

follows from the N = 2 case of Lemma C.2. Hence we arrive at ( 15. 23[) . □ 

Thus far we have assumed y is real. When we fix x G M 3 , however, it is easy to see 
from the defining formula (15. 6p and Prop. 4.4, combined with the bounds (I5.2p - (I5.3I) . that 
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J 3 is well defined for y G C 3 whose imaginary parts are sufficiently close. Indeed, we easily 
obtain an estimate 

\M x iV)\ < ■ / t t wp m / dzexpi-a^ + z 2 )\m(y l + y 2 - 2y 3 )/2) 

sin ( Trim (yi -y 2 )/Reb) J R2 

x sinh(q(2:i - z 2 )) sinh(7lm (y x - y 2 ){z l - z 2 )/2) ^ 
Uj=i nLi cosh( 7 (z fc - X j)) 

and the integral occurring here converges when |Im (yi — y 2 )\ and |Im (yi + y 2 — 2y 3 )\ are 
smaller than Re b. This also entails that J 3 is holomorphic in y in the latter domain. 

In fact, it turns out that a more refined bound gives rise to a larger y-domain, allowing 
also a subdomain of D% for the dependence on (b,x). This hinges on Lemma C.3, which 
enables us to calculate the integral occurring in a suitable majorization of Js(6;x,?/). To 
arrive at the latter, we need to start from a quite different representation for J 3 , which is 
only well defined on a subset of D 3 . The details now follow. 

To begin with, we introduce the notation 

1 3 1 3 

x 3 = -^2xj, Y 3 = -^2yj, Xj = Xj-X^ yj = yj -Y 3 , j = 1,2,3, (5.29) 
i=i i=i 

and recall that the representation (15. lip is well defined for any (b,x) G D 3 . By contrast, 
we need to restrict (6, x) to 

D r 3 = {(b,x) G S a x C 3 I |Im%| < a-Re&/2, j = 1, 2, 3} C D 3 , (5.30) 

for the formula 

e iay 3 (x 1 +x 2 +x 3 ) r 

J 3 (b;x,y) = / dzl 3 (b;x,y,z), y G M 3 , (5.31) 

1 J(M.+iImX 3 ) 2 

to be well defined and valid. This representation can be rewritten as 

Jz{b;x,y) = ey^{?>iaX 3 Y 3 )J r 3 {b]X]y), (5.32) 

j;{b;x,y) = l dsW 2 (b;s)Sl(b;x,s)J 2 (b;s,( yi -y 3 ,y 2 -y 3 )), (b,x) G D r 3 , y G R 3 , 

(5.33) 

where we used ( 14.16]) . It is to be noted that J' 3 depends on x and y only via the differences 
xi - x 2 , x 2 - x 3 and y x - y 2 , y 2 - y 3 . 
Next, we introduce the domain 

V 3 = {(b, x, y) G D r 3 x C 3 I |Im ( Vj -y k )\< Re 6, j, k = 1, 2, 3}, (5.34) 

and the maximum function 

u(x) = max |Imx,-|. (5.35) 

i=l,2,3 J 

Now we are prepared for the analog of Prop. 4.4. 
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Proposition 5.4. The function Js(b;x,y) is holomorphic in (b,x,y) on T> 3 . For any 

(b, x, y) e V 3 with 

Re (xj - x k ) ^ 0, Im (yj - y k ) 0, 1 < j < k < 3, (5.36) 

we have 

\J 3 (b;x,y)\ < C(b, /i(x)) exp(-3alm (X 3 F 3 )) 

T,res 3 (-) T exp ( - a ^=1 Re (^0)) Im %) 
rii< m <„<3 sin(vrlm (y n - y m )/Re b) sinh(7 Re (x m - x n )) ' 

where C is continuous on S a x [0, a — Re 6/2). 

Proof. In view of (15. 32 p . we need only estimate J 3 . From Prop. 4.4 and the bounds (15. 2ft - 
(15. 3p we readily deduce 

C(b /i(x)) f 

\ J l{ x ,y)\ < • 7 t 7 ' w-p ,x / dsexp(-a(si + s 2 )Im(yi + y 2 -2y 3 )/2) 

sm(7rlm (?/! - y 2 )/Re&) y R2 

sinh(a(si - s 2 )) sinh(7im (?/i - y 2 )(si - s 2 )/2) 



n,=i IL=i cosh( 7 (s fc - Rexj)) 



(5.38) 



We proceed to rewrite this in such a way that Lemma C.3 may be invoked. To this end 

we set 

S 2 = ( Sl + s 2 )/2, sj^Sj-Sz, j = 1,2. (5.39) 
Consider now the auxiliary integral 

a /, \ f r sinh(7(si - s 2 )) / „ v-^/ X 

A 2 (t,u)= ds 3 2 — exp ( 7 5 2 2^(M3-Mi)J 

^(-) T exp (7^(1)^2 (5.40) 



TSS 2 

Clearly, v4 2 ( 7 Rex, 2Im?//Re&) is equal to the integral in the bound (I5.38P times a factor 
—2. The integrand is symmetric under the interchange s\ s 2 , and the sum is antisym- 
metric. Thus we may replace the sum by twice the summand with r = id. Taking next 
s — > z/j, this yields 

2 ^ ^ sinh(;zi - 2 2 ) 
r3 tt2 



A 2 (t, u) = — j dz- , 

7 Jr2 n j= i IL=i cosh(z fc - t j/ 

2 

x exp (^Z 2 ^2(u 3 - UjYj exp ((zi - Z 2 )(u 2 - m)). (5.41) 



j'=i 

Comparing to C 2 (u,t) (1C3 1 1) . we deduce 



A 2 (t,ii) = -^C 2 (tt,t). (5.42) 

7 2 



23 



Therefore, the bound (15.381) can be rewritten as 



\j;(x,y)\ < - , ^Y {x) \ m M C 2 (2Imy/Re6,7Rex). (5.43) 
7 Z sm(7rlm (7/2 — 2/i) / He 0) 

We can now appeal to Lemma C.3 to deduce the holomorphy assertion. Also, assum- 
ing (I5.36p . we can combine ( 1C.34j) with 



2 3 

^2 Re ( £ r(j))lm (y 3 - Vj) = Re (^(i)) Im (V3 - %) 
i=i j=i 
3 3 

= - Re (x r(j) ) Im ^ = - Re (x r(j) )lm y jt (5.44) 

j=i i=i 

to obtain (Q7|l . □ 

Thus far we have obtained holomorphy properties of J3 by employing two contours 
equal to the same horizontal line. We conclude this section by allowing more general 
contours, so as to obtain a larger holomorphy domain. However, we still keep both 
contours equal and shall also restrict attention to real y. Indeed, as will transpire shortly, 
even with these restrictions there are novel complications arising for N = 3 (as compared 
to N — 2), and we shall not strive for optimal results. In fact, we do not try and obtain 
a counterpart of Prop. 4.6, since even with such a result available, it would be a major 
undertaking to extend the N = 4 case beyond Props. 6.1-6.4. 

Just as for J 2 , we can deform the contours to show that J 3 remains holomorphic for 
(b, x) in a larger domain than D 3 . We restrict attention to deformations of the contours for 
which both contours remain equal to the same contour C, defined by functions f(s) from 
the space C w delineated above ( IB. 181) . The choice of contour is such that it separates 



the upward and downward pole sequences ( 15. 5p . Also, to stay clear of the W/Vpoles 
for |Im6| > 0, we need / to be such that when \s\ increases by |Im 6| , then \f{s)\ increases 
by less than |Re6|, cf. (1B.18j) . This can be achieved in two distinct ways. We can either 



require 

\f{s)\ < \Reb\/2, VsGM, (5.45) 
or we can restrict the slope of C: 

\f{s)\ < |Refe|/|Im6|, Vs e R. (5.46) 

Next, the contour C should be further restricted to ensure that the J 2 -factor in the 
integrand (15.11) remains holomorphic when Z\ and z 2 vary over C. For b G (0, 2a) this can 
be readily achieved, since in that case J 2 is holomorphic for Kezi 7^ Re z 2 . For Imb ^ 0, 
however, we need |Im(zi — z 2 )\ < 2a — Re b whenever Re (z\ — z 2 ) equals ±Im&. Again, 
we can ensure this by requiring either 

\j{s)\ < |2a-Re6|/2, Vs G R, (5.47) 

or 

\f(s)\ < |2a-Re6|/|Im6|, Vs G R. (5.48) 
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Accordingly, in the following analog of Prop. 4.5 we obtain an extended holomorphy 
domain for J3 by allowing only contours for which either 

\f{s)\ < m(Reo)/2, Vs G R, (5.49) 

or 

|f (s)| <m(Re6)/|Im6|, Vs G R, (5.50) 

where m(d) is given by 

m(d) = min(2a — d, d), d G (0, 2a). (5.51) 
Specifically, we define two domains 

Di = {(6,x) G S a x C 3 | |Re(x m -x n )| ^ |Im6|, 

max |Im(xj — Xfc)| < 2a — Re6 + m(Re6)/2}, (5.52) 



ee {(6,ar) G 5 a x C 3 | |Re(x m -x n )| > |Im6| > 0, 

|Im (x m — x n ) Im b\ < (|Re (x m — x n )\ — |Im&|)m(Re 6)}, (5.53) 

where 1 < m < n < 3 and k = 1, 2, 3. 

Proposition 5.5. Let y G R 3 . TTien t/ie function J 3 (b;x,y) is holomorphic in (b,x) on 
the domain 

£>g = D 3 U £>£ U D$. (5.54) 

Proof. First assume (6, x) G D3. We start from the representation (15. lip . Then the first 
pole in the upward/downward sequences (15. 5p has an imaginary part that is below/above 
the contour R + i<p(x) by a distance less than m(Re6)/2, so we need only use down- 
ward/upward triangular indentations whose heights are smaller than m(Re&)/2. (Note 
that we can make disjoint indentations, since |Re(x m — x n )\ 7^ |Im6|.) This ensures 
that Z\ — Z2 stays away from the HVpoles and the J2-argument stays in D 2 - Hence 
holomorphy in results. 

Next, let (b, x) G Df. The restriction on the real parts implies that the numbers 
Rc distinct. By permutation symmetry, we need only define the 

piecewise linear contour C for the case Re £3 < Re £2 < Rexi. Then we have 

Rex 3 + |Im6|/2 < Rex 2 - |Im6|/2, Rex 2 + |Im&|/2 < Rexi - |Im&|/2. (5.55) 

First, we connect the pairs of points 

P± = x m ± |Im6|/2, m = l,2,3, (5.56) 

with (horizontal) line segments. Then we connect P 3 + , P 2 ~ and P 2 + , P{ with line segments. 
The restrictions on D% imply that the slopes of the latter satisfy (I5.50p . Finally, we 
connect P 3 ~ to (-R, 0) and P x + to (R, 0), where R > is chosen sufficiently large so that 
the slopes s± of the two line segments also satisfy \s±\ < m(Re6)/|Im6|. Completing C 
with the intervals (—00,—R) and (R, 00), we stay in the region where the integrand is 
holomorphic. Hence holomorphy of J 3 in follows as well. □ 
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6 The inductive step 



In this section we use the validity of Props. 5.1-5.4 as the starting point for the inductive 
step N — 1 — > N. More specifically, Props. 6.1-6.4 below reduce to the former for N = 3, 
and our induction assumption is that they hold true if we replace N by N — 1. This 
has the advantage of expository conciseness, but the minor drawback that we need to 
invoke Prop. 6.3 with iV — > N — 1 to obtain Prop. 6.1. As it turns out, our account in 
Section 5 can be followed to a large extent. Accordingly, we sometimes skip details when 
their general- N version will be clear from Section 5. 

The integrand arising when constructing J N from Jn-i is given by (cf. (I2.16P ) 

I N (b; x, y, z) = W N ^{b] z)S N (b; x, z)J N ^(b; z, (y x -y N ,..., y N ^ - y N )). (6.1) 

From flB.6j) we obtain the bound 

N N-l r (h _ , 

i^v)i<nn ^ B ;;;?| i)r ™ 

as long as (6, z k — Xj) G V r for j — 1, . . . , N and k = 1, . . . , N — 1. Also, from (1B.13P we 
have 

|Wjv_i(6;«)| < ] [ C w (b) sink 2 (^(zj - z k )), Imz x = • • • = Imzjv-i- (6.3) 

i<j<k<N-i 

In order to estimate the factor Jn-i in ( 16. ip . we assume that the imaginary parts 
of zi, . . . , Zn-i are equal to some c G R, and replace N by N — 1 in Prop. 6.3 below. Then 
the difference function cIn-i(z) (given by (16.201) and (16. 8ft ) vanishes and (pN-i(z) (16.7ft 
reduces to c. Hence we obtain the bound 

| J N -i{b) z, {y x -y N ,..., y N _ x - y N ))\ 

N-l 



<C(b)exp[-acy2(yj-y N )] 1 [ . " — rr, (6.4) 

» 1 - LJ - smh 7 (z m — z n )) 

<m<n<N-l \ l\ m nj ) 



where C is continuous on S a . Together with the bounds (16.2ft and ( 16. 3ft . this implies 
that In has exponential decay when one or more of the quantities |Rezi|, . . . , |RezAr_i| 
diverge, uniformly for c varying over bounded intervals. 

We note that the poles of the kernel function S N (b;x,z) are located at (cf. (lA.lip - 
(|A~12]) ) 

Zk = Xj±(ib/2-ia — Zt m ), k = 1, . . . , N - 1, j = l,...,N, l,meN. (6.5) 

In addition, by the induction assumption, we may make use of Prop. 6.1 after taking 
N N — 1, and hence conclude that Jat_i(5; 2, (yi — ?/jv, . . . , yN-i — Vn)) is holomorphic 
in Dn-i- Requiring at first x G R , it now follows that the integrand Jjy is regular for 
2 G R^ 1 , so that 

e iay N (xi^ \-x N ) /" 

■M&l y) = / rfz/7v(&; 2/, 2), (6, x, i/jG^xfx (6.6) 
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is well defined. 

It is clear that Jj^(b]X,y) (as initially defined by ( 16. 6p ) extends to a function that is 
holomorphic in 6 and Xj for 6 G S a and |ImXj| < a — Re 6/2, j = 1, . . . , N. Also, the 
above bounds imply that we are allowed to shift all contours R up and down by the same 
amount, provided none of the poles (16.51) are met. In this way we can recursively extend 
the holomorphy domain. In order to make this precise, we set 

<Mx) = Im {x h + x jN )/2, x G C N , (6.7) 

where the indices ji and ]n are such that 

lmx jN < Imx jjv _ 1 < • • • < Imi i2 < Imx^, {j 1 , . . .,j N } = {1, . . -,N}. (6.8) 

Then we have the following general- N version of Prop. 5.1. 

Proposition 6.1. Let y G M. N . Then the function Jjv(6; x, y) is holomorphic in 

D N = {(b,x) e S a x C N \ max llm ( Xj - x k )\ < 2a - Reb}. (6.9) 

l<j<k<N J 

Proof. Let us fix (6, x) G D^. Then we have 

Imxj + a — Re 6/2 > <p N (x), Imxj — a + Re 6/2 < cf> N (x), j — l,...,N. (6.10) 

It follows that the contour M. + i<pN{x) stays below/above the upward/downward G-pole 
sequences (I6.5p . Thus we can continue Jn to the x- value in question by simultaneously 
shifting all contours, without passing any of these poles. In this way we obtain the 
representation 

e iay N (x!-\ \-x N ) r 

J N (b;x,y) = — — — — — / dzI N (b;x,y,z), (b,x) G D N , y G R N , (6.11) 

[1\ - L)l jQSL+idNix))*- 1 

whence holomorphy in is plain. □ 

In our next proposition we prove that Jn is a joint eigenfunction of the 2N A AOs (ll.2ip . 
As in Prop. 5.2, we need to shift some of the N — 1 contours in the proof. To show that 
this does not cause problems with the contour tails, we can no longer rely on the bounds 
(16. 3p and (16 .4p . since they only hold for Imzi = ■ ■ ■ = Imz^^i. 

However, it follows from the bounds (IB. 14[) ( TB. 15[) that the rate at which Wn^i(z) 
diverges as |Re(z m — z n )\ — > oo, 1 < m < n < N — 1, is independent of Im (z m — z n ). 
Moreover, by the induction assumption, we may invoke Prop. 6.3 for N — > N — 1. As long 
as |Im(z m — z n )\ < 2a — Re 6, this entails that on closed subintervals the decay rate as 
|Re (z m — z n )\ — > oo in the bound on the JAr_i-factor in the integrand In ( 16. ip does not 
depend on Im (z m — z n ) either. Combined with the bound ( 16 .2p on the kernel function S N , 
these arguments show that the contour shifts are legitimate. 

Proposition 6.2. Fixing y G and letting (b,x) vary over 

D N ,s = {(b,x)EDN\(b,x-ir ] )EDN,yv^[0,ai} N , Re6G(0,a s )}, (6.12) 

we have the eigenvalue equations 

A${x)J N {x, y) = Sl N) (e 5 (2 Vl ), . . . , e s (2y N ))J N (x, y), (6.13) 

where k = 1, . . . , N, and 5 = +, — . 
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Proof. The six steps in the proof of Prop. 5.2 can be readily adapted, so we only indicate 
the changes. The first two steps yield 



e s {2ky 



e iayjv(^iH \-x N ) 



N 



(N-l)\ 

x I dz x ---\ dz N _iW N _ 1 (z)A ( k N s \x)S s N (x,z)J N ^i(z, (y x -y N ,...,y N ^ - y N )), 

(6.14) 



C_ Jc 



as the generalization of (I5.14p . 

For the third step, we invoke the kernel identity (12. 8p and the symmetric function 
recurrence (11.281) to deduce that we need only show that the integrals 

dzi-- • / dz N -iW N -i(z)J N -i(z, (yi -y N ,-- - , 2/jv-i - Un))A ( ' 1 ^ 1 \-z)s1 1 (x, z), 

'C- JC- 

(6.15) 

are equal to 

S i N ' 1] M 2 (l/i - Vn)), e 5 {2(y N . l - y N ))) 

x (N — 1)! exp(— iay N (%i + h x N ))J N (x, y), (6.16) 

for 1 = 1,...,N- 1. 

Fourth, we fix I and use ( 11.211) to write (16.151) as 



dz x ■ • ■ / dz N -iW N -. 1 (z)J N - 1 (z, {y x -y N ,.. . ,y N -i ~ Vn)) 

C_ JO- 

x II S5{Z 1 " + f } II exp(ia_,^j4(x, z). (6.17) 

IC{1,...,N— 1} mSJ mGY 
|/|=i 

Changing variables z m — > z m — ia_s, m G /, we rewrite (I6.17P as 

ssi.z m - ia-g - z n + ib) 



\i\=i "4i 

x II ex P(-* a -<5d*jV^v-i(20^iV-i(2, (Vi - Vn,-- ■ ,Vn-i - Vn)), (6.18) 

and then use (11.201) to see this equals 



j2 (n / dz ™\ (n / dz n )w N ^(z)s 

_,...,N-l}\meI Jc -+ ia -S J \n0 Jc - J 



n(x, z) 



/C{1,...,7V-1} 
1/1=/ 



x Yl ~^~7 — ~ — T"^ II ea P(~ 4 ' a -^*»)«^-i(2, (z/i - Vn,-- • ,2/at-i - Vn))- (6.19) 

nil 
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We claim that when we simultaneously shift the contours C_ up by a^s for all variables 
z n , n I, we do not meet poles. (The paragraph preceding the proposition implies that 
these shifts do not create problems at the tail ends.) Assuming this claim is valid, we note 
that the integrands are regular for Zi, . . . , Zjv-i G C_ + ias, an d since the shifts ensure 
that the N — 1 contours are equal in all integrals, we may use the eigenvalue equation 
(I6.13P with iV — > N — 1 and k — > I (the fifth step). In the sixth step, we shift the contours 
back to M, hence arriving at ( 16.161) . 

It remains to prove the claim. Obviously, we stay in the holomorphy domain of S^ N 
while shifting the contours. Using Prop. 6.1 with N — > N — 1, we see this is also true 
for Jn-i- Now we fix / C {1, . . ., N — 1} such that |/| = I and first choose a_s = <V 
When shifting the contours for z n , n ^ /, we only meet the simple poles of the factors 
1/G(z m — z n + ia — ib) with m G / and n ^ I at z m — z n = ib, but they are matched by 
zeros of the factors Ss(z m — z n — ib). Moreover, the simple poles arising from ss(z m — z n ) 
at z m — z n = and z m — z n = ia>i (encountered when a s = a/) are matched by zeros of 
the factors G{z m — z n + id) in Wn-i(z). It follows that the integrands are regular for 
z m G C- + ia s , m G J, and z n in the strip \m.z n G [—ai/2, —ai/2 + a s ], n ^ J, as claimed. 

Choosing next a^s = a-u the simple poles of 1/G(z m — z n + ia — ib), m G I,n ^ /, 
located at z m — z n = ib + ika s for ka s < ai — b, are matched by zeros of ss(z m — z n — ib), 
and the simple poles resulting from ss(z m — z n ) at z m — z n = ika s for ka s < a>i are matched 
by zeros of G{z m — z n + ia). Hence the integrands are also regular for z m G C_ + ia h 
m G /, and ln\z n G [— a;/2, ai/2], n ^ /, as claimed. □ 

The generalization of Prop. 5.3 involves the function <Pn{x) 116. 7ft and the distance 
function 

cLn{x) = Im (xj 1 — Xj N ), x G C , (6.20) 
where ji and jn are given by (16. 8p . 

Proposition 6.3. Let y G WL N . For any (b,x) G we have 



|Jjv(&;x,y)| < C(b, d]^(x))exp I— a 



N N-l 

<p N (x)^y,j +yN^2 ( Im x i* ~ M x )) 

J=l k=2 

FT Re(x m -x n ) 
i< m <t<iv smh ^ Re ^- x ™)) 



where C is continuous on S a x [0, 2a — Re b). 

Proof. From the representation (I6.1ip we have the majorization 

exp ( -ay N J^f=i Im^) 
|Jjv(x,y)|< V 



(AT-1)! 



x / dz\W N -i(z)S 9 N (x,z)J N -i(z,(yi-y N ,...,y N -i-y N ))\. (6.22) 

By the induction assumption, (I6.2ip holds true with N replaced by N — 1. Applying this 
to the factor Jn-i of the integrand, we observe that d^^i^z) = and (f>N^i(z) = (p^{x). 
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Using also the bounds (16. 2 j) and (I6.3p . we now infer 

(AT AT-1 
-a <f) N {x) y^Vj + Vn^2 ( Im x ^ _ 

x/ ds n - o smh( 7 ( gm - Sn) ) n n Cr(6; *; t ^ ~ * j) . (6.23) 

.V- 1 ,^_iix 7 , cosh(7(s fc -Rex,)) 



\<m<n<N-l j- 



Just as in the proof of Prop. 5.3 (cf. (I5.26P ). we see that the C r -product is bounded above 
by a function C(b, cIn(x)) that is continuous on S a x [0, 2a — Re 6). Using Lemma C.2 to 
compute the remaining integral, (I6.2ip results. □ 

As in the N = 3 case, we can allow complex y, provided we restrict attention to a 
subdomain of D N for the dependence on (b,x). First, we introduce general- N notation 
that deviates slightly from the notation (15.291) for the case N = 3: For a given vector 
v E C M , M > 1, we define V M G C and v {M) G C M by 

1 M 

Vm =mJ2 v ^ ^ M) = Vj -V M , J = 1,...,M. (6.24) 

3=1 

Then we need to work with 

D r N = {(b,x) eS a xC N \ \lmx<j N) \ < a -Re 6/2, j = 1, . . . , N} C D N , (6.25) 
for the representation 

p iay N (x\-\ hx N ) 



J N (x,y)= — / dzI N (x,y,z), y e R JV , (6.26) 

to hold true. 

We now use (I6.26j) as a starting point for concluding that Jat(x, y) equals the product 
of a center-of-mass factor exp(NiaX N Y N ) and a function J r N {x, y) depending only on the 
differences Xj — Xj + i, yj — yj+i, j = 1, . . . , N — 1. The resulting representation is important 
in its own right, and it will enable us to control the y-dependence. Specifically, we claim 
that we have 

J N (x, y) = exp(NiaX N Y N )J r N (x, y), (6.27) 
J^{x,y) = / dsW N -i(s)S i N (x {N) ,s)J N _ 1 (s,(y 1 -y N ,...,y N „ 1 -y N )), 

(iV - i)\ J R N-1 

(6.28) 

where (6, x) e D r N and y e M. N . 

As before, our proof of this claim proceeds by induction on N, noting (16. 27j) — (I6.28j) 
reduce to (I5.32l) - (l5.33p for N = 3. First, we observe that the integral on the rhs of (16.261) 
can be rewritten as 

rfsPUAr_i(s)5^(x (Ar) , s)J N -i((si + X N , . . . , sat_i + X N ), (y 1 -y N ,-- ■ ,Vn-i - Vn))- 

(6.29) 
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Next, we take N — > N — 1 in (I6.27p -f l6.28p . to deduce from the induction assumption 
Jn-i{{si + X N , sat_i + X N ), (y x -y N ,..., y N _ x - y N )) 

= exp I iaX N 2j(% - y N ) I Jn-i(s, (yi -yN,-- -,Vn-i ~ Vn))- (6.30) 

Using the identity 

N N-1 

y N Xj + X N (Vj ~ Vn) = NX N Y N} (6.31) 

3=1 3=1 

the claim now readily follows. 
Setting 

V N = {(&, x, y) e D r N x C N \ \lm ( Vj - y k ) \ < Re b, j, k = 1, . . . , N}, (6.32) 

/ijv(x) = max |Ima;,- |, (6.33) 

j=l,...,N 

we are in the position to state and prove the arbitrary- iV version of Prop. 5.4. 



Proposition 6.4. The function J^{b;x,y) is holomorphic in (b,x,y) on V N . For any 
(6, x, y) G V N with 

Re( Xj -x k )^0, lm( yj -y k )^0, 1 < j < k < N, (6.34) 

we have 

\J N (b;x,y)\ < C(b,/j, N (x))exp(-Nakn(X N Y N )) 

E (Te5jv (-rexp(-aRe(xW) -a(lm^))) 



rii< m <„<7v sin(7rlm (y n - y m )/Reb) sinh(7Re (x m - x n )) 
where C is continuous on S a x [0, a — Re 6/2). 



(6.35) 



Proof. It is clear from (I6.27P that we only have to estimate J r N . The induction assumption 
amounts to (I6.35P holding true when we replace N by N — 1. Combining the resulting 
bound on Jn-i with ( 16. 2 p and (I6.3j) . we obtain via a straightforward computation (using 
the general- iV counterpart of f )5.44p ) 

|JM«,»)I< C(b ' Mx)) 



( N - l V-Ui< m <n<N-i sm (^ lm {y n -y m )/Reb) 

f , ril< J <fc<^-l Sinh ^( S i ~ S k)) ( V^t / 

Jr n - 1 , , . , coshf'YiSfc — Rear- ')) \ 



1 11^=1 Ilfc=i cosh(7(s fc -Rex} ')) \ j=i 

/ JV-2 

^ (-) T exp aj] 4oi 1)lm (^-i ~~ ) ■ ( 6 - 36 ) 



x 

rG5]v_i \ i=l 
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Next, we introduce 



A U \- f A IIl< 3 -<fc<JV-l Slnh (7(3j -Sk)) ( ^ ( 

A N -x{t,v) = / ds - N — — — exp jS N -i > (u N -Uj, 

Jrn-i [7 11 cosh( 7 s fc - tj) \ p[ 



=S N -! \ j=l J 



X 

tGSjv_i 

It is plain that A^-i^Rex^ , 2lmy /Reb) equals the integral in the bound (16. 36 p . Now 
we fix r G Sjv-i and consider the corresponding summand on the rhs of (I6.37p . Changing 
variables Sj — > s T -i^ and using antisymmetry of the sinh-product, we obtain the sum- 
mand with r = id. Hence we may replace the sum by (N — 1)! times this summand. We 
simplify the resulting product of exponential functions by using the identity 

N-l N-2 JV-1 

S N -i £Otv - uj) + £ sJ^Viv-i - Uj) = £ Sj(u N - uj), (6.38) 

j=l j=\ j=i 



which is readily verified. Taking s — > z/7 and comparing the result to (1C.31|) . we deduce 

(N - IV 

A N ^(t,u) = v / Cjv-i(M,t). (6.39) 
The upshot is that the bound ( 16.36!) becomes 

\J r A^y)\ < ^TTT C ^^ff C N ^(21my/Reb,jRex^). 

(6.40) 

Using again the general- N version of (15.441) . the assertions now follow from Lemma C.3. 

□ 

We point out that the bound (16.351) is symmetric in y, a feature we also expect 
for Jtf(x,y), cf. the next section. Another point of interest is its specialization to real y. 
Using (p22jl -( 1023]l we see that this is of the form 

\J N (b;x,y)\<C(b,fi N (x))exp(-NaY N lmX N ) TT . ^ ^ m - %n) y G R N . 

i<J£<n smh( 7 Re(o; m -x n )) 

(6.41) 

It should be noted that for iV > 2 this does not coincide with the bound ( I6.2ip obtained 
in Prop. 6.3. Moreover, the latter holds true on Dn, whereas we need (b,x) G D r N to 
arrive at (16.411) . 

We conclude this section by pointing out that since \J^(b; x,y)\ is manifestly non- 
negative, it follows from (I6.35P that the function on the rhs is positive. This corollary 
of Prop. 6.4 is by no means obvious, even for N = 3. Rephrased in terms of the func- 
tion Qjsi-iiu^t) given by (1C33j) . it implies 

Qjv-i(M) > 0, \uj-u k \<2, l<j<k<N, u,teR N . (6.42) 
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7 Outlook 



In this paper we have focused on the first steps in the recursive scheme for construct- 
ing joint eigenfunctions Jn of the AAOs A k ^ (given by (12.161) and (ll.2ip . resp.). More 
specifically, we have concentrated on establishing holomorphy domains and uniform decay 
bounds that were sufficient for proving that the scheme does provide well-defined func- 
tions Jn that satisfy the expected joint eigenvalue equations. There are numerous aspects 
that remain to be investigated. Below we briefly discuss a few that we plan to study in 
future articles, associated with conjectured features of the joint eigenfunctions. In the 
free cases b = a± these features all follow from the results in Section 3, whereas for N = 2 
they can be gleaned from }Ruillj . 

Sn x Sn invariance: We expect that Jat(x,?/) is invariant under permutations of 
both the variables x = (xi, . . . ,xn) and the variables y = (yi, . . . , 2/jv)- We add that 
the former symmetry property is immediate from the defining formula (I2.16p . whereas 
the latter is far from clear. We also recall that the bound (I6.35P exhibits the expected 
symmetry. 

Self-duality: In the free cases b — a$, 5 — ±, it is plain from Theorem 3.1 that upon 
multiplication of Jn^s] y) by the factor rii<j<fc<7v s s(Uj ~ Vk) / S-sillj ~ Vk) we obtain a 
function that is symmetric under the interchange of the variables x and y. More generally, 
for the function 

J N {b;x,y) = J N (b;x,y) ] [ G{8(yj - y h ) + ia - ib), (7.1) 

±<j<k<N 5=+,- 

we conjecture the self-duality property 

J N (b;x,y) = J N (b;y,x). (7.2) 

Parameter symmetries: Since the hyperbolic gamma function G(a + ,a_; z) is in- 
variant under the interchange a + -H- a_, it is clear from the defining formula (12. 16ft that 
JN{a+,CL-,b;x,y) has this parameter symmetry, too. Moreover, as we noted in the in- 
troduction, the Hamiltonians H^ s are invariant under the substitution b — » 2a — b. We 
conjecture that their joint eigenfunctions 

x, y) = G{ib - ia) N ^ N -^ 2 W N (b; x) l l 2 J N (b; x, y)W N (b; yf/\ (7.3) 

are also invariant under this substitution: 

7"jv(2a - b;x,y) = F N {b;x,y). (7.4) 

Global meromorphy: Although we have only considered holomorphy properties in 
this paper, we do expect that the recursive scheme yields globally meromorphic joint 
eigenfunctions. Moreover, the free cases and the N = 2 case suggest both the location 
of the poles and upper bounds on their orders. Specifically, we expect that the product 
function 

V N (b;x,y) = J N (b;x,y) ] J ] [ E(5(x j -x k )+ib-ia)E(5(y j -y k )+ib-ia), (7.5) 

l<j<k<N 6=+,- 

has a jointly analytic extension to all (b,x,y) G S a x C N x C N . (The .E-function is an 
entire function related to the hyperbolic gamma function by G(z) = E(z)/E(—z), from 
which the location and order of the E- zeros can be read off, cf. [Rui99j .) 
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Soliton scattering: It is a long-standing conjecture that the particles in the relativis- 
tic Calogero-Moser systems of hyperbolic type exhibit soliton scattering (conservation of 
momenta and factorization). In more detail, provided b e (0,2a), we expect 

JF N {b;x,y) ~ C N ^ W ( 

o-G5iv i<j<k<N 

a-\j)«r-\k) 

x n (~ u ( y j 

l<j<k<N 
a- 1 (j)>a- 1 (*) 

for Xj — Xj + i — » oo, j = 1, . . . , N — 1, with CV a constant and the scattering function u{z) 
given by (11.171) . 

Orthogonality and completeness: We conjecture that for all b G [0,2a], the func- 
tion A/jvJ-jv(o; x, y) (with A/jv a normalization constant) yields the kernel of a unitary 
operator on L 2 (F^,dx), where is the configuration space 

F N = {— oo < Xjv < • ■ ■ < Xi < (7.7) 

In the free cases this unitary operator amounts to the antisymmetric version of Fourier 
transformation, cf. Section 3. 



-u{Vj -Vk)Y /2 

- Vk)y 1/2 exp(zax ■ <r(y)), (7.6) 



Appendix A. The hyperbolic gamma function 

In this appendix we collect properties of the hyperbolic gamma function G(a + , a_; z) we 
have occasion to use. Their proofs can be found in Subsection III A of |Rui97], but for 
the G-asymptotics we need the sharper bounds obtained in Theorem A.l of |Rui 99j . We 
fix 

a+, a_ > 0, (A.l) 

and suppress the dependence of G on a + , a_ whenever this will not cause ambiguities. 

The function G(z) can be defined as the unique minimal solution of one of the two 
analytic difference equations 

G(, + W2) 
G[z — ias/2) 

that has modulus 1 for real z and satisfies G(0) = 1 (cf. (II .9p for the notation used here); 
it is not obvious, but true that this entails that the other one is then satisfied as well. It 
is meromorphic in z, and for z in the strip 

S = {z e C | |Im(z)| < a}, a=(a+ + a_)/2, (A.3) 

no poles and zeros occur. Hence we have 

G(z) = exp(ig(z)), z e S, (A.4) 

with g(z) holomorphic in S. Explicitly, g(z) has the integral representation 

( \ f°° d y ( sin2?/2 z \ 

g(a +J a„;z)= — , — — , z e S. A.5 

Jo y \2sinh(a + y)sinh(a_y) a+a^y J 
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From this, the following properties of the hyperbolic gamma function are immediate: 

G(—z) = 1/G(z), (reflection equation), (A. 6) 

G(a_, a + ; z) = G(a + , a_; z), (modular invariance), (A. 7) 

G(\a + , Aa_; Xz) = G(a + ,a_; z), A G (0, oo), (scale invariance), (A. 8) 



G(a + , a_; z) = G(a + , a_; — z). (A. 9) 

Defining 

z kl = ika + + ila_, k, I G N = {0, 1, 2, . . .}, (A.10) 
the hyperbolic gamma function has its poles at 

z = z^ t = —ia — z k i, k,leN, (G— poles), (A. 11) 

and its zeros at 

z = z£ t , z kl = ia + Zki, fc,l€N, (G— zeros). (A. 12) 
The pole at — ia is simple and has residue 

lim (z + ia)G(z) = — (a+a_) 1/2 . (A.13) 

2— >— ia 2,71 

More generally, for a given (fc , ^o) G N 2 , the multiplicities of the pole zj^ an d zero z^ q1q 
are equal to the number of distinct pairs (k,l) G N 2 such that z^ = . In particular, 
for a + /a_ ^ Q all poles and zeros are simple. 

Finally, we specify the asymptotic behavior of G(z) for Re (z) — > ±oo. To this end we 
introduce a function f(a + , a_; z) by setting 

G(a+, a_; z) = exp ( =)= i (y + a-2 2 /4 + f(a+, a_; z)) ) , ±Rez > a^, (A. 14) 

where 

* = £ (- + -)■ < A15 > 

24 \a_ a+ / 
It then follows from Theorem A.l in |Rui99] that we have 

\f(a + , a_; z)\ < C(a + , a_, Imz) exp(— aa s \Rez\/2), \Rez\>ai, (A. 16) 

where C is continuous on (0, oo) 2 x K. Here and in the main text, we find it convenient 
to formulate bounds that are uniform on compact subsets of a given set S in terms of 
positive continuous functions on S, generically denoted by C. 
It follows from the above that we have 

\G(z)\ < C G (z)exp(almz\Rez\/2), VzeP G , (A. 17) 

where Cq is continuous on the domain 

V G = {z G C | z i -i[a, oo)}, (A.18) 

and satisfies 

C G (z) = 1 + 0(exp(-aa s |Rez|/2)), \Rez\ -»■ oo. (A.19) 

Here, the implied constant can be chosen uniformly for Imz varying over a bounded 
interval. 
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Appendix B. Bounds on G{z—ib /2) / G(z+ib /2) and w(b] z) 

In the main text we have occasion to use a complex parameter h whose real part is 
restricted to (0, 2a). This entails that the downward pole sequence of the function G(z — 
ib/2) is located in the lower half z-plane and the upward pole sequence of 1/G(z + ib/2) 
in the upper half z-plane. With b varying over the strip 

S a = {z £ C | Re z £ (0, 2a)}, (B.l) 

the ratio function 

r(6; z) ee G(z - ib/2) /G(z + ib/2), (B.2) 

plays a crucial role. It is expedient to estimate it in two different ways, depending on 
reality properties. Below and in Sections 4-6, we often use a new parameter 

7 = a Re6/2 = ^^. (B.3) 
a + a_ 

Lemma B.l. Letting b £ (0,2a), we have 

\r{b;x)\ < C r (6)/cosh(a6x/2), Vs £ R, (B.4) 
with C r (b) a continuous function on (0,2a). Moreover, defining 

V r ee {(b,z) eS a xC\±z- ib/2 <£ -i[a, oo)}, (B.5) 

we have 

\r(b;z)\ < C r (M/cosh(7Rez), V(6, z) £ V r , (B.6) 
where C r (b, z) is a continuous function on T> r satisfying 

C r (b,z) = l/2 + 0(exp(-aa s |Rez|/2)), |Rez| -»■ oo, (B.7) 

wrf/i the implied constant uniform for (6, Im z) varying over compact subsets of S a x R. 

Proof. From ( 1A.14I) we obtain 

r(b;z) =exp ( =f [q;6^/2 + if(z - ib/2) - if (z + ib/2)]), ±Rez > \lmb\ + a t . (B.8) 



In view of ( 1A.16I) this entails 

2 cosh (7 Re z)r (6; z) = 1 + 0(exp(-aa s |Re ^|/2)), |Re^|->-oo, (B.9) 

uniformly for (6, liaz) in compacts of S a x R. Hence, fixing b £ (0,2a), we readily 
obtain ( IB. 41) . Since r(b; z) is analytic in T> r , the second assertion easily follows from (IB. 91) . 

too. □ 

Next, we consider the weight function 

w(b;z)= Yl G(5z + ia)/G(8z + ia-ib). (B.10) 
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By flA.9j) it is real for b and z real, and positive for z G R*. Furthermore, it is even and 
holomorphic in 



V w = {(b, z) G S a x C | ±z + zfo ^ -i[0, oo), ±2 ^ -z[2a, oo)}, 
and on account of ( 1A.6I) and ( 1A.13I) it satisfies 

4vr 2 



w(b; z) 



G(ib-iayz 2 + 0(z 4 ), b G S a , z -»■ 0. 



(B.ll) 



(B.12) 



Just as for r(6; z), we obtain two different bounds on the weight function, cf. Lemma B.l. 
Lemma B.2. Letting b G S a , we have 

\w(b;x)\ <C w (b)smh 2 (-fx), \/x G R, (B.13) 
mi/i C w (b) continuous on S a . Moreover, we have 

\w{b; z)\<C w (b,z) cosh 2 (7 Re z), V(b,z)eV w , (B.14) 
where C w (b,z) is a continuous function on T> w (IB.llj) . which satisfies 

C w (b, z) = 4 + 0(exp(-aa s |Re z|/2)), |Re z| -> 00, (B.15) 
uniformly for (b,lmz) varying over compacts of S a x R. 
Proof. Letting ±Re z > \lmb\ + a + a;, we get from (IA.14I) and (IA.6[) 



w(b; z) = exp ( ± abz — i ^[/(^ + + /(^ + <5(z6 — ia))] 

«=+,- 



(B.16) 



□ 



Hence the assertions readily follow by using (1A.16j) and (IB. 121) . 

In Section 5 (below Prop. 5.4) we employ contours C in the z-plane of the form 

C = { s + ic + if(s) I s G R, c G R, / G £}, (B.17) 

where £ is the space of real-valued, continuous, piecewise linear functions with compact 
support. There are further restrictions on the functions ensuring that the contours stay 
away from singularities. In particular, focussing on w(b] Z\ — z 2 ) for b G S a and Z\,z 2 G C, 
we need ±(zi — Z2) + ^ to stay away from the pole sequences in —i[Q, 00). Therefore, we 
restrict attention to the subset C w of functions in C that satisfy 



± ( Sl -s 2 ) = lmb^ ±(/( Sl ) - /(s 2 )) £ (-00, -Re6], b G 



(B.18) 



Note that the pole sequences arising for ±(z\ — z 2 ) G — i[2a, 00) are always avoided. 
Indeed, s\ = s 2 implies z\ = z 2 , since Z\ and z 2 belong to the same contour. 
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Appendix C. Three explicit integrals 

The integral in the first lemma enables us to render the scheme explicit in the free cases b = 
a±, cf. Section 3. The next two lemmas may be viewed as corollaries. Lemma C.2 is a 
key tool to bound the joint eigenf unctions recursively for real y, yielding holomorphy 
in the variable x in suitable domains, cf. Props. 6.1 and 6.3. Lemma C.3 enables us to 
include holomorphy properties in the variable y, cf. Props. 5.4 and 6.4. 

Lemma C.l. The integral 

I N (p,t)= dzexp(tz -p) ~ — -, (C.l) 

Jn N \[j=l l h=i cos Htj ~ z k) 

converges absolutely for all p e C N and t e <C N+1 satisfying 

|Impi|<2, i = l,...,N, llmtj-l < 7r/2, j = 1, . . . , N + 1. (C.2) 

For such p, t, with in addition pi ^ ; i = 1, . . . , N, and tj ^ tk for l<j<k<N + l, 
it is given by 

N 1 N 

^ fet) ^^^ l J w+Ii ^ T E(-^exp( ^ E^). (C3) 

Proof. The first assertion readily follows from the bound 

Y[ swh( Zj -z k )\ <exp((JV-l)(|«i| + -.- + |^|)), WzeR N . (C.4) 

l<j<k<N 

In order to prove the explicit evaluation formula (1C.3I) . we first observe that by Fubini's 
theorem the value of Jjv is independent of the order of the multiple integration. Next, 
starting with the zjv-integral, we note the equality 



sinh(7rp Ar /2) / dz N exp(iz N p N 



n^L/ sinh^ - z N ) 
Ilf=i cosh(tj - z N ) 



1 f llyLi 1 sinh(zj — z N ) 
= o 1^ s dz N exp(i(z N -isn/2)p N )— WTl — -. (C.5) 

2 s =+,- U j= i cosh(t i - z N ) 

Changing variable z^ — > z^ + isn/2, the right-hand side becomes 

1 f nf=i cos ^( z j — z n) 
~o/^ s dz N exp(iz N p N ) : -. (C.6) 

1 S T^_ JR-isn/2 H j=1 smh(tj - Z N ) 

In the strip |Imzjv| < tt/2, the only poles of the integrand are the simple poles located at 
Zjsr = tj, j = 1, . . . , N + 1. Invoking Cauchy's theorem, we arrive at the sum 

N+l TtJV-1 Uf j- \ 

• f J H cosh( Zj -t lN ) 

m > exp (it lN p N — — — — . C.7 

^ n^smh^-^) 
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With the ZAr-integral now determined, we obtain 



N+l 

ITC 

'TV 



m ST exp(iti N p N ) t(1n) 



sinh(7rpjv/2) ^ FL^i* sinh(t^ - t,J" 



where 



4-1 = / dz i ■ ■ ■ dz N-i exp ( i z jPj J pj^^Sfcl rT-^ k ~ • (C9) 

jk 1 *- 1 \ y y n fc =i coBh(tj - z k ) 

Since ijv^i is of the same form as ijv-i> we can proceed recursively to compute the 
remaining integrals. In this way we end up with 



±± sinh(7r Pj /2) ^ P(h,...,l N ;t) 1 1 ' ' 



JV 

Z7T 

'TV " ' 



where 



P(h,...,l N ;t)= Yl sinhfew-^)"' II sinh^-tjj, (C.ll) 

and the prime signifies that the sum extends over 

= 1, . . . , JV + 1, ^l N ,...,h^l N ,...,l 2 . (C.12) 

To each such set of indices l\, . . . ,1^, we now associate a uniquely determined permutation 
r G tSjv+i by requiring lj = r(j), j = 1, . . . , N. Clearly, we have 

P(T(l),...,T(N);t) = s(r) J] sinh^-4), (C.13) 

i<i<jfc<Ar+i 

where s(r) equals 1 or —1. The product on the right-hand side of (1C.13|) yields a function 
that is antisymmetric in t, whereas I^(p,t) is manifestly symmetric in t. Thus s(r) must 
be equal to (— ) r up to an overall sign. Letting r = id, we see this sign equals (— ) N , so 
that we arrive at (lC.3j) . □ 

The absolute convergence of ( 1C.1I) for all (p, t) G C w x C N+1 in the polystrips ( 1C.2I) 
implies holomorphy of 7(p, i) in the resulting product domain, but there is no useful closed 
formula for the coincidence limits of the right-hand side of (1C.3|) . Also, we can freely take 
Pj-partials under the integral sign in view of the exponential decay of the integrand for 
Zi — > ±oo. In particular, in Section 5 we shall encounter an integral 

B*{t)= [ dz jg ~f Sinh ^ ~ Z2 \ = lim(-id Pl + id P2 )I 2 (p, t), (C.14) 

whose direct calculation from ( 1C.3j) would already be unwieldy. 
In order to avoid this chore, we introduce 

F(p,t) = J2(-) Te Mir{t)-pl (C.15) 
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where both p and t belong to C 3 . By antisymmetry of F in p and t, its power series 
expansion is of the form 

F{ P ,t) = G{ P ,t) Y[ (p m -Pn){tm-t n ), (C.16) 

l<m<n<3 

where G(p, t) is symmetric in p and t, and invariant under interchange of p and t. Fur- 
thermore, we have 

3 

G(p, t)=c +J2 c jkPjtk + 0((p m t n ) 2 ), (C.17) 
j,k=i 

where the notation will be clear from context. It is not hard to verify that Co equals —i/2, 
so that we can now deduce 



F((p 1 ,p 2 ,0),t) = -'-p 1 p 2 (p 1 -p 2 ) 11 (t m -t n )(l + R(p u p 2 ,t)), (C.18) 

Km<n<3 



where the remainder R vanishes for pi = p 2 = 0. When we now use the representa- 
tion (IC.lgp to rewrite I 2 {p,t) as specified by (1C.3I) . then a moment's thought shows that 
the integral (IC.14p is given by 

The following lemma details the arbitrary- N version of this integral, which we need 
in Section 6. 



Lemma C.2. Let t E C N+1 , with |Imt,-| < ir/2, j = 1, . . . , N + 1. Then the integral 

ril<m<n<7V ( Zm ~ Z n) sfnll(2: m — Z r , 

nK 1 nf=i cosh(t i - Zk 

trn tr. 



B N {t) = j dz ~ ^ N ~ +1 N — , (C.20) 



is given by 



B N {t)=2 N • , V ( C - 21 ) 



Proof. We shall proceed in a slightly different manner compared to the N = 2 case 
considered above. However, the starting point is the same, namely, the function 

F N (p,z) = {-YeMir(z)-p), (C.22) 

t€SW 

where p, z E C N . Reasoning as before, we find that its power series expansion is of the 
form 

f N \ 

c ikPj z k + 0({p m z n ) 2 ) J JJ (p m - p n )(z m - 2 n ). (C.23) 

j,fc=l / Km<n<JV 
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From this we deduce 



F N ((pi,.. . ,ptv_i,0),z) 

= VnPi ■ ■ -Pn-i Yl & ~ II ( Zm ~ Zn "> ' i 1 + R & z ">) ' ( C - 24 ) 



l<j<fc<iV-l l<m<n<N 



where the remainder R(p, z) vanishes at p = 0. To determine the constant t]n, we note 
that the monomial {piZ\) N ~ l (p2Z2) N ~ 2 ■ ■ ■ Pn-\Zn-i only occurs in the expansion of 

1 (i PlZl + ■■■ + iPN z N ) N ^-^ 



(N(N -l)/2) 

N(N-l)/2 

= i N(N-l)/2 y- {PlZl) n ■ ■ ■ {PNZ N y N , Q25 , 

Ji,— ,Jjv=0 
ji + -+j N =N(N-l)/2 

Consequently, we have 

VN = (N-mN-2)\..-2V (C ' 26) 
Now from flC.23j) we see that the integral B N (t) is given by 

B N (t) = hmn N — -. (C.27) 

Ul<m<n<N(Pm ~ Pn) 

On the other hand, requiring at first tj ^ tk for 1 < j < k < N + 1, and using r(£) • a(p) = 
(a -1 o r)(t) ■ p, o G Sn, we deduce from (1C.3I) that the numerator can be written 



A? 



-Z7T 

sinhf 

It follows from fl(J.24[) that the rhs of this equality is of the form 



1 Km<n<Jv+l 



(-2<)^w n (Pi-f*) n sinbTt _n ' (i+ ^ (p ' o) ' (c - 29) 

l<j<k<N l<m<n<N+l ^ m n ' 

where again the remainder R(p,t) vanishes at p — 0. Substituting this expression for 
the numerator and observing that (— i) N\t]n+i = we arrive at ( 1C.21j) . Since both 
the integral and the rhs of ( 1C.21I) are well defined when some coordinates coincide, the 
evaluation formula follows for all t in the specified polystrip. □ 

To state the final lemma (which we invoke for Props. 5.4 and 6.4), we introduce 

1 - 

Zn =nY. z ^ zeRN - ( c - 3 °) 
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Lemma C.3. The integral 

f rii<i<fc<A r sinh^ — Zk) /v-^ \ 

C N (u, t) = / dz ~ — exp f 2^ Zj{u N+l - Uj) ) , (C.31) 

J* N U j= i 1L=i cosh(t 3 - - z ft ) V =1 / 

converges absolutely for all (u, t) in the domain 

C N = {u,te C N+1 | \Re{u N+l -Ui)\ < 2, i = 1,...,N, |Imt,-| < vr/2, j = 1, . . . , N + 1}. 

(C.32) 

TTie quotient function 

Q N (u,t) = C N (u,t)/ [ J (uj-u k ), (C.33) 

l<j<k<N 



is holomorphic in Cn- For (u, t) G Cat, wift m addition Uj ^ u^+i for i = 1, . . . , N, and 
tj tk f or l<j<k<N + l, we have 

= n n < - — ^ n • ur , s 

" sm(7r(wjv+i - Uj)/2) ± ^^ N+1 smh(i j - * fc ) 

iV 

x S (-) T exp(^t T(i )K +1 - M ,)). (C.34) 

tgSjv+i j=l 

Proof. Comparing ( 1C.31I) to ( 1C.1I) . we deduce 

C N (u,t) = I N (p(u),t), p(u)j = i{uj - u N+ i), j = l,...,N. (C.35) 

Hence the first assertion and the explicit evaluation ( 1C.34j) follow from Lemma C.l. The 
evaluation reveals that C^{u,t) is antisymmetric under permutations of ui, . . . ,Un, so 
that holomorphy of Qn in Cn readily follows. □ 

At the end of Section 6 we are in the position to deduce a remarkable positivity feature 
of QN(u,i) from the use of Lemma C.3 for obtaining the bound (16.35p . cf. fl6.42j) . Under 
the stronger assumption Uj G (un+i — 2,un+i), j = 1, ■ ■ ■ ,N, this feature can also be 
derived directly from Lemma C.3 (using finite induction). However, due to the special role 
of itjv+i, a direct proof of (16.421) in the setting of this appendix would be quite laborious. 
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